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ABSTRACT 



The averaged probability of detecting a planetary companion of a lensing star during a 
microlensing event toward the Galactic center when the planet-star mass ratio q = 0.001 
is shown to have a maximum exceeding 10% at an orbit semimajor axis near 1.5 AU 
for a uniform distribution of impact parameters. This peak value is somewhat lower 
than the maximum of 17% obtained by Gould and Loeb (1992), but it is raised to 
more than 20% for a distribution of source-lens impact parameters that is determined 
by the efficiency of event detection. Although these probabilities, based on a signal to 
noise (S/N) detection criterion, are model and assumption dependent, the fact that they 
change in predictable ways as functions of the orbit semimajor axes but remain robust 
for plausible variations of all the relevant galactic parameters implies that they are 
representative of real values. In addition, the averaging procedures are carefully defined, 
and they determinine the dependence of the detection probabilities on several properties 
of the Galaxy. The probabilities for other planet-star mass ratios can be estimated 
from an approximate scaling of ^fq. A planet is assumed detectable if the perturbation 
of the single lens light curve exceeds 2/ (S/N) sometime during the event, where it is 
understood that at least 20 consecutive photometric points during the perturbation are 
necessary to confirm the detection. S/N is the instantaneous value for the amplified 
source. Two meter telescopes with 60 second integrations in I-band with high time 
resolution photometry throughout the duration of an ongoing event are assumed. The 
probabilities are derived as a function of semimajor axis a of the planetary orbit, where 
the peak probability occurs where a is approximately the mean Einstein ring radius of 
the distribution of lenses along the line of sight. The probabilities remain significant 
for 0.6 < a < 10 AU. Dependence of the detection probabilities on the lens mass 
function, luminosity function of the source stars as modified by extinction, distribution 
of source-lens impact parameters, and the line of sight to the source are also determined, 
and the probabilities are averaged over the distribution of the projected position of the 
planet onto the lens plane, over the lens mass function, over the distribution of impact 
parameters, over the distribution of lens and sources along the line of sight and over the 
/-band luminosity function of the sources adjusted for the source distance and extinction. 
The probability for a particular impact parameter and particular source / magnitude 
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but averaged over remaining degenerate parameters also follows from the analysis. In 
this latter case, the extraction of the probabilility as a function of a for a particular q 
from the empirical data from a particular event is indicated. 
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1. Introduction 

There are several searches for extrasolar plan- 
ets using radial velocity techniques that together 
have discovered more than 50 planets with minimum 
masses less than 13 Jupiter masses (Mj) (e.g., Marcy, 
Cochran & Mayor 2000) . Astrometric searches so far 
have been unsuccessful, but new astrometric satellites 
(Full-sky Astrometric Mapping Explorer (FAME), 
Space Interferometry Mission (SIM), Global Astro- 
metric Interferometer for Astrophysics (GAIA)) with 
precisions approaching micro-arcseconds (Lattanzi et 
al. 2000) and the Keck Interferometer (Swain et al. 
2000), which may do almost as well, should soon add 
additional planets. Transit searches have yet to dis- 
cover a new planet, but a planet discovered by radial 
velocity techniques has been successfully observed in 
transit (Charbonneau et al. 2000; Henry et al. 2000). 
Continuing transit searches should lead to more dis- 
coveries. All of these search techniques except as- 
trometry favor the discovery of close, massive plan- 
ets, and all are restricted to planetary periods that 
are only slightly greater than the time the program 
has been in operation or less. The one astrometric 
program with a long time base (Gatewood, 1991) but 
no detected planets has too few target stars to be 
definitive. The program initiated by G. Marcy, while 
having a long time base, has reached the necessary 
precision of 3 m/sec to detect distant Jupiters and 
Saturns only a few years ago, but the current pro- 
gram has a sufficient number of stars to constrain the 
statistics of distant Jupiters similar to the one in our 
own solar system in about 5 years (G. Marcy, D. Fis- 
cher, private communications, 2000). 

Mao and Paczyhski (1991) pointed out that a plan- 
etary companion of a star acting as a gravitational 
lens amplifying the light of a more distant star (mi- 
crolensing) could perturb the otherwise smooth light 
curve in an easily observable way and thereby reveal 
its presence. There are many examples in the lit- 
erature of microlensing light curves perturbed by a 
planet along with descriptions of the single lens am- 
plification (e.g. Peale, 1997). The information about 
the planet from the perturbed light curve alone is 
usually limited to the planet-star mass ratio q and 
the projected separation of the planet from the star 
in terms of the Einstein ring radius, 
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where G is the gravitational constant, M is the mass 
of the lensing star, c is the velocity of light, D OL and 
D os are the observer-lens and observer-source dis- 
tances respectively, z — D l/D$, and £ = D os /Ds 
with Z?8 = 8 kpc being the distance to the center 
of the Galaxy. Statistically, the projected star-planet 
separation is a measure of the semimajor axis of the 
planetary orbit (Peale, 1997). 

Part of the search for microlensing events has 
been directed toward the Galactic bulge, where the 
high concentration of stars in the telescope field al- 
lowed a reasonable rate of event detection and where 
more than 300 events have, in fact, been detected 
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princeton.edu/6glc/ ; http://www.lal.in2p3.fr/rccherche 
/eros/erosa.html). Gould and Loeb (1992) showed 
that the probability of detecting a planet with q = 
0.001 during a microlensing event approached 17% if 
the planet's projected separation from the star was 
comparable to the Einstein ring radius R E . A planet 
was assumed detectable if the perturbation exceeded 
5% of the instantaneous flux of the unperturbed single 
lens light curve. The probability was averaged over a 
uniform distribution of lenses along the line of sight 
(LOS) to a source at the Galactic center. The aver- 
aged position of the lens was thus halfway to a source 
at the galactic center, where R E » 4.1 AU for a solar 
mass star meant that a Jupiter would have nearly the 
maximum probability of detection. 

This high probability of detecting a planetary com- 
panion during a microlensing event generated an en- 
thusiastic advocacy for establishing an intensive, ded- 
icated microlensing search program from the ground 
(Beichman et al., 1996). Initially more modest mi- 
crolensing search programs (PLANET (Probing Lens- 
ing Anomalies NET work), Albrow et al. 1998, and 
MPS (Microlensing Planet Search), Rhie et al. 2000) 
have progressed to different degrees of maturity and 
are beginning to place meaningful constraints on the 
occurrence of Jupiter mass planets around M dwarfs 
(Albrow et al. 2000). The enthusiasm was reinforced 
by the special status noted for very high magnification 
events, where the probability of detecting a planet is 
unity for a range of projected planet-star separations 
in the lensing zone surrounding the Einstein ring ra- 
dius (Griest and Safizadeh 1998). The advantage of a 
microlensing search over the other search techniques 
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is that the planets in long period orbits could be dis- 
covered without watching for a decade, and statistics 
would start accumulating immediately. The disad- 
vantage, of course, is that the lenses are so distant 
that there could be no followup studies of the planet 
after the event. But the knowledge of the frequency of 
occurrence of planetary systems and the need to know 
these statistics rapidly for planning future planetary 
search missions are sufficient to keep microlensing as 
a most important search technique. 

This work develops a probability of planet detec- 
tion during microlensing events that is based on the 
photometric signal-to-noise ratio (S/N) for the am- 
plified source, with the goal of determining the de- 
pendence of this probability on plausible variations 
in Galactic parameters while establishing robust av- 
eraging techiques over the distributions of Galactic 
parameters. The probability is determined as a func- 
tion of the planet orbit semimajor axis for a planet- 
star mass ratio of q = 0.001, and an approximate 
scaling law for other mass ratios is indicated. From 
some points of view, it might be more desirable to ex- 
press the probabilities for planetary masses instead of 
the planet-star mass ratio. However, the microlens- 
ing calculations are naturally expressed in terms of 
q, and any attempt to extract a probability in terms 
of the planetary masses would have increased already 
lengthy calculations by an order of magnitude. As it 
is debatable whether mass or mass ratio is the more 
interpretable quantity, the small return is not worth 
the considerable extra effort. To determine a depen- 
dence of the detection probabilities on Galactic pa- 
rameters, we must necessarily average the probabili- 
ties over distributions of quantities, such as lens mass, 
that remain inherently unknown and over other quan- 
tities such as source-lens impact parameter that are 
not a priori known. Details of the dependence of the 
detection probability on most of the parameters in- 
volved are shown as these averages are calculated. In 
this process we show that the S/N criterion leads to 
maximum averaged probabilities of detection that are 
comparable to or even greater than that obtained by 
Gould and Loeb for plausible assumptions of Galac- 
tic parameters, event geometries and observing tele- 
scopes. Perhaps more importantly, these probabilities 
remain robust as the parameters are varied within rea- 
sonable ranges. We should emphasize to the reader 
that the probability is not that of detecting a planet 
during any random event, but only the probability 
of detection if the lens has a planetary companion. 



Hence, the success of any search program depends on 
the frequency of planetary occurrence, but null results 
will be meaningful in the sense that if the planets are 
there, they will be detected with reasonable probabil- 
ity. 

An S/N criterion has been used previously to de- 
termine empirical planet detection probabilities for an 
individual event (Gaudi and Sackett, 2000; Albrow, et 
al. 2000a; Albrow et al. 2000b). A best fit single lens 
light curve for a point source-point lens model is de- 
termined from a given event data set by minimizing 
Xpspl with free parameters t n , t E} u min , F , f B . The 
time of minimum separation between source and lens 
is to; the event time scale t E = R E /v±, where v± is 
the relative transverse velocity of source and lens pro- 
jected onto the lens plane; the source-lens impact pa- 
rameter is u mm = jmin/jE, with 7 min being the mini- 
mum angular separation of source and lens during the 
event and j E = R E /D OL ; the unamplificd flux density 
from the source is F ; and the ratio of the blended flux 
density from unresolved stars to the source flux den- 
sity is f B . To check if a binary lens gives a better fit to 
the data, one assumes a fixed planet-lens mass ratio 
q, a fixed projected planet separation b = r/R E , and 
a projected source trajectory through the lens plane 
inclined at angle a relative to the lens-planet line and 
again minimizes \ 2 with the same set of free parame- 
ters. A planet is assumed detected if Ax 2 (q 7 x p ,a) = 
xlspL-X 2 {q,x P ,ot) > A X 2 thresh = 100, where the lat- 
ter value is considered sufficiently large to be a secure 
better fit to the data (Gaudi and Sackett, 2000). To 
determine the efficiency of detecting a, q,b planet for 
a particular event, the above process is repeated for a 
uniform distribution of < a < 2tt, and light curves 
which would have A\ 2 — X 2 (?, x p , a) — x 2 PS pl greater 
than some threshold would indicate that the planet 
would have been detected. The detection efficiency 
(probability) 

1 f 27T 

e{x p , q) = — Q[A X 2 {x p , q, a) - A X 2 thresh }da 

"""" (2) 
is the fraction of trajectories for which x 2 ( x p> q, a) ex- 
ceeds the threshold value for detection for this (x p , q), 
where <d(x) is a step function. The probability in 
terms of the orbit semimajor axis in AU and planet 
mass m is estimated by assuming representative val- 
ues of the lens mass and the lens and source distances. 

This procedure confirms that planet detection dur- 
ing a microlensing event has a significant probabil- 
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ity for individual events. We show in Section 7 how 
e(x p , q) can be averaged over the degenerate param- 
eters in the Einstein ring radius weighted by the dis- 
tribution of lenses and sources along the LOS and 
over the mass function (MF) yielding a probability of 
detection P'(a,q) applicable to that particular event 
with its determined source magnitude and normalized 
impact parameter. The Gaudi-Sackett procedure can- 
not be used to obtain the overall averages sought here 
for arbitrary events. Here we have no data sets on 
which to establish a x 2 criterion, so the reference light 
curve will be that of the single lens without the planet, 
and the perturbations by the presence of a planet will 
be relative to this single lens curve. Since in the real 
world one has only best fit single lens light curves 
and best fit binary lens light curves for comparison, 
both of which tend to suppress any perturbations, our 
probabilities will be somewhat over estimated. This 
overestimate will be at least partially compensated 
by our procedure which tends to underestimate the 
probabilities. 

We justify below the assumption that a planet will 
be detectable if the fractional perturbation of the light 
curve and hence the mean of a set of points during 
the perturbation exceeds 2/(S/N) for the magnified 
source. Implicit in this assumption is that the pertur- 
bation be observed sufficiently long that enough pho- 
tometric points can be obtained to establish the path 
of the photometric deviation with more than 99% con- 
fidence. In this sense, the S/N for the entire set of 
anomalous points will be VJJ greater than that for an 
individual point, where A/" is the number of photomet- 
ric points during the perturbation. We argue below 
that Af = 20 will establish a detection with greater 
than 99% confidence if the true mean of the set of 
perturbed points is displaced at least 2/ (S/N) from 
the unperturbed single lens light curve. Our neglect 
of the non-zero angular size of the source will limit 
q > 10~ 4 , where perturbation durations (2^/qt E ) will 
be a significant fraction of a day. So for almost all 
binary geometries where the planetary perturbation 
will exceed 2 /(S/N), the duration of the perturbation 
will be sufficient to establish a reliable mean of the 
light curve during the perturbation — especially since 
the sampling frequency will most likely be drastically 
increased during the perturbation. 

The S/N used as the detection criterion is de- 
veloped in Section ||, where the S/N from photon 
statistics and sky noise is reduced by a constant fac- 
tor to account for systematic and other unmodeled 



noise and thereby approximate two independent em- 
pirical determinations of S/N for real systems. As we 
shall need a luminosity function (LF) to determine 
the fraction of the sources that are visible at any dis- 
tance D os , we also introduce in this section the LF 
in J-band along with its adjustment for the source 
distance and extinction thereto. We have chosen the 
Holtzman et al. mass function (MF) for the lenses 
but the Zocalli et al. (2000) LF for the sources. Al- 
though Holtzman et al. (1998) determine an /-band 
LF for the Galactic bulge on the way to deriving their 
MF, it spans only about 9 magnitudes. But Zocalli 
et al. (2000) have obtained a deeper bulge LF in 
J-band for the lower main sequence and have supple- 
mented their main sequence LF with the bulge giant 
LF obtained from Tiede et al. (1995) and Frogel and 
Whitford (1987) to obtain a more complete LF for 
the bulge that spans 15 magnitudes. The straight line 
segment representation of the main sequence part of 
the Holtzman et al. /-band LF has the same slope as 
the LF given in Section || within the uncertainties, so 
we keep the MF of Holtzman et al. (1998) that was 
derived from their LF, since programs incorporating it 
had been written earlier and since the MFs of several 
other studies also have an index of 1 in the M star re- 
gion (e.g. Basu and Rana, 1992; Gould, Bahcall and 
Flynn, 1997 if unresolved binaries are accounted for 
(Gould, private communication, 1998)), but use the 
more complete /-band LF developed from the J-band 
LF of Zocalli et al. (2000) in Appendix A. Zocalli et 
al. (2000) derive a slightly steeper MF in the M star 
region than Holtzman et al. (1998) (indcx=1.3 vs 1). 
The change of the index from 1 to 1.3 in the M star 
region has only a small effect on the detection proba- 
bility. We choose to illustrate the effect of increasing 
the MF index in the M star region by the relatively 
large index of 2.2. 

We introduce the fundamental idea of the proba- 
bility of detection based on an S/N criterion in Sec- 
tion ||. Here contours of the projected position of 
the planet onto the lens plane, within which the frac- 
tional perturbation of the single lens light curve ex- 
ceeds 2 /(S/N) in magnitude, are shown to sweep out 
an area in the lens plane as the source passes by the 
lens at an impact parameter u m i n = r m i n /R E . The 
extent of the the area grows with increasing q, with 
increasing S/N , and with decreasing u m i n . We shall 
assume q = 0.001 for most of the calculations, but 
we also show a completely averaged probability as a 
function of the planet semimajor axis for q = 0.0001. 
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(The S/N determined from Section |] is a function of 
the /-band magnitude of the amplified source for a 2 
meter telescope and 60 second integration times.) If 
the projected position of the planet is anywhere inside 
this area, the planet will cause a fractional perturba- 
tion of the single lens light curve exceeding 2 /(S/N) 
for a time of order 2^Jqt E sometime during the event 
and be detected (e.g. Gould and Loeb, 1992). The 
probability of detection is just the probability that 
the projected position of the planet is within the area 
swept out by the contours. (The amplification of the 
source for various configurations of the lens-planet bi- 
nary system and the generation of contours of con- 
stant perturbation of the single lens light curve used 
extensively in this section are developed in Appendix 
B.) Also illustrated in this section is the reason for 
the unit probability of detection in the annular "lens- 
ing zone" region straddling R E for high amplification 
events (Griest and Safizadeh, 1998).) 

The probability of the projected planet position is 
developed in Section |] and Appendix C. To average 
over unknown orbital inclination and planet orbital 
longitude, a planet is assumed to be uniformly dis- 
tributed over a sphere of radius a, and the probabil- 
ity F(r) that the projected separation of the planet 
from the lensing star is between r and a is deter- 
mined therefrom (Gould and Loeb, 1992). F(r) is ex- 
pressed in terms of x p — r/R E and a/R E , as this is the 
normalization under which the microlensing amplifi- 
cation is expressed. Then F(r) — ► F(x pi a, (, z, M). 
The explicit display of M, z and C in F(x p , a, C, z, M) 
follows from R E (Eq. (|l])), and the degeneracy of 
these variables in any microlensing event makes nec- 
essary an average over z from observer to source 
weighted by the distribution of lenses along the LOS 
and an average over M weighted by the MF to 
yield F(x p ,a,() (Appendix C). The average over £ 
weighted by the distribution of sources along the LOS 
must follow an average of the probability over the 
LF. F(x p , a, C) is the probability of finding the pro- 
jected position of the planet between x p = r/R E and 
a/R E averaged over the lens MF and the distribu- 
tion of lenses along the LOS from the observer to 
the source, and —dF(x p ,a,Q)/dx p is the probability 
density. The Zhao (1996) bulge model is combined 
with the Bahcall-Soneira (1980) disk model to de- 
termine the spatial distribution of lenses. We have 
already pointed out our choice of the Holtzman et 
al (1998) MF above. The probability of detection is 
now initially developed for fixed a and £ for a par- 



ticular source /-band apparent magnitude. Eleven 
values 0.6 < a < 10.0 are carried through the remain- 
ing steps to obtain a final averaged probability as a 
function of a. 

In Section || for q ~ 0.001 and given I-band mag- 
nitude of the source, a and u m i n , we determine the 
probability P(a,C,I,u m i n ) of detecting a planet of 
semimajor axis a, when the normalized source dis- 
tance is C) the source magnitude is / and the source- 
lens impact parameter is u m i n . This is accomplished 
by determining the boundaries of the area swept out 
in the lens plane by the contours of the 2/ (S/N) frac- 
tional perturbation of the single-lens light curve and 
integrating the probability density of the planet pro- 
jection over the area, where the latter simply uses 
differences in F(x p , a, £) for two pairs of values of x p 
defining the boundaries of the area at stepped inter- 
vals along the source trajectory. P(a, £, I,u m i n ) is 
determined for a distribution of impact parameters 
< Umin < 1, and the results averaged over u m i n 
to yield P(a,(,I), which is the probability of detec- 
tion for the given q = 0.001, semimajor axis, source 
distance and source /-band magnitude averaged over 
the distribution of planetary orbital inclinations and 
phase, over the distribution of lenses along the LOS, 
over the MF, and over the distribution of u m i n . These 
probabilities are determined for a grid of values of 
C, /. It is also shown here how the peak of the detec- 
tion probability moves toward smaller values of a as 
the source moves closer to the observer, which reflects 
its following the mean value of R E . 

For each £ the span of likely values of / is deter- 
mined according to the LF as modified by the dis- 
tance to and extinction of the source with / < 21 
assumed to be the largest magnitude for which use- 
ful high precision, high time resolution photometry 
could be carried out. [For data sets for which image 
differencing (e.g., Alcock et al. 2000) or other tech- 
niques have not been employed to pick up many of 
the events with fainter sources, the number of usable 
faint sources will be reduced from that implied by the 
real LF with / < 21. Reducing the fraction of sources 
with low luminosities would actually increase the av- 
eraged probabilities, since the brighter stars will have 
a higher S/N. (See Fig. 0.)] The P(a,(,I) is then 
averaged over I min < I < 21 weighted by the LF for 
each £ to obtain P(a, £) for each £. I m i n corresponds 
to the brightest star kept in the LF; it decreases as 
the source distance C is decreased or the extinction is 
reduced. Illustrations of the dependence of the detec- 
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tion probability on source distance and source magni- 
tude for averages over the LF at each source distance 
are also given in Section ||. The dependence of / on £ 
is why we delayed averaging over the latter until the 
average over the LF could be completed. 

For 0.1 < £ < 1.2, P(a, £) is averaged over the dis- 
tribution of visible sources along the line of sight in 
Section || to yield P(a), where "visible" means I < 21 
as discussed above. P(a) is the probability of detect- 
ing a planetary companion of a lens with q = 0.001 
and semimajor axis a (for 11 values of a) during a mi- 
crolensing event averaged over all of the Galactic and 
geometric parameters. The same galactic model used 
for the distribution of lenses is used for the distribu- 
tion of sources. The lower limit on £ > is chosen to 
avoid a singularity in the expressions, where a source 
at £ < 0.1 is extremely unlikely in any case. The up- 
per limit is chosen 1.2 since extinction and blending 
will make the number of sources that are usably visi- 
ble negligibly small at larger £. It is shown here that 
for a LOS toward Baade's window, the peak probabil- 
ity of detection exceeds 10% for a near the mean value 
of R E of almost 2 AU, if u min is distributed uniformly, 
but it exceeds 20% for a distribution determined by 
the observational efficiency for the events, which has 
a higher proportion of small u min . 

To clarify the series of averages, we write an equa- 
tion for P(a) that includes all of the averaging op- 
erations except the initial one over the planetary or- 
bital inclination and phase, where the latter is accom- 
plished simply by the assumption of uniform distribu- 
tion over a sphere of radius a. 



P(a) 
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(3) 



where @[f(x p , a, £, z, M)} is a step function whose ar- 



gument > [given in Eq. (CI)] is a condition for the 
reality of F(x p , a, £, z, M); n' L (z, M) is the fraction of 
lenses in volume z 2 dzAtt and in mass range dM of 
those along the LOS to the source with Ail being a 
representative solid angle in the field; S{u l min , I) is the 
area in the lens plane corresponding to the particular 



u min an d source 7-band magnitude within which a 
planet would cause a perturbation of the light curve 
exceeding 2/(S/N). The double integral on the far 
right is F(x p ,a, £) so —dF(x p , a,0/dx p is the prob- 
ability density, which is integrated over the area in 
the next integral to the left to yield P(a, £, I, u min)- 
The sum of P(a, £, J, u % min ) over the k values of u m i n 
divided by k gives the average over the impact param- 
eters to yield P(a, £, I). LF(I, £, A) is the fraction of 
sources per unit I magnitude adjusted for distance 
and extinction A. This integral over the LF yields 
P(a, £). Finally, n' g (£) is the fraction of visible sources 
of those along the line of sight, and the integral over 
C yields P(a). 

Section ^ shows how the detection probability 
changes as various parameters are changed. It de- 
creases as the extinction is increased or if the detec- 
tion criterion is increased from 2/(S/N) to 3/(S/N). 
Changing the LOS changes the distribution of the 
probability over the values of a by changing the mean 
R E . Here we show how increasing the fraction of Mi- 
type stars in the MF drives the peak in the probability 
to smaller values of a, while decreasing the probabil- 
ity for large a and increasing the probability for small 
a. We have already pointed out above the conver- 
sion of the model independent e(x p ,q) to the model 
dependent P'(a, q) at the end of Section 7. 

The properties of the averaged planetary detection 
probability are summarized in Section ^ where it is 
also shown how the detection probability, averaged 
over everything except the LF, varies with the source 
magnitude. Also shown here is the detection proba- 
bility as a function of the planetary semimajor axis 
for a particular event, where the source magnitude 
and the impact parameter are assumed known. A 
discussion follows in Section |9|, where the robustness 
of relatively large detection probabilities over a wide 
range of planetary semimajor axes to variations in 
the various parameters is emphasized, while the de- 
pendencies on assumptions and other caveats and ne- 
glected processes are pointed out. The unique aspects 
of a microlensing search for planets and the robust- 
ness of the probabilities of detection for all plausible 
variations in Galactic parameters are additional jus- 
tifications for the continued pursuit of ground based 
microlensing searches for planets in spite of the expo- 
nentially growing success of radial velocity searches 
and proposed space based microlensing searches. 
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2. Signal to noise ratio. 

The part of S/N appropriate to photometry of a 
star that depends on just photon statistics and sky 
noise can be written (up to the factor 1/4 in paren- 
theses) 



S 
N 



r y/A T Et 



1 



y/A a $, kv + (j) sta 

314 x 10-° AI DVt 
y/2.23 x 10- 8 + 10-°- 4/ 



(4) 



where (j> 3tar (photons/cm 2 sec) is the photon flux den- 
sity from the source star, A T — nD 2 /A (cm 2 ) is 
the telescope aperture area (D = aperture diame- 
ter), E is the overall efficiency of the system, § sky 
(mag/arcsec 2 — > photons/(cm 2 sec arcsec 2 )) is the sur- 
face brightness of the sky, A s — 7rs 2 /4 (arcsec 2 ) is 
the seeing disk (s = seeing defined as the FWHM 
of the point spread function) and t (sec) is the in- 
tegration time. From the lines just preceding Eq. 
(A7), the /-band flux density from Vega is 1.56 x 
10~ 6 ergcur 2 sec -1 leading to <f> ttar = 6.28x 10 5 10"°- 4/ 
photons cm _2 scc _1 for magnitude I with a square 
profile for the /-band filter 0.13/im wide at 0.8/xm 
(Ivega = 0). The clear night sky brightness at Cerro 
Tololo Interamerican Observatory in /-band ranges 
from 19.2 to 19.9 magnitudes arcsec -2 (Walker, 1987) 
from full to new moon. We choose an average sky 
brightness of 19.6 magnitudes arcsec -2 , s = 1.4 arc- 
sec and E = 0.2 to yield the final form of Eq. (0) 
without the factor (1/4). 

Eq. (ji|) without the factor 1 /4 does not yield values 
of S/N close to those empirically determined because 
of instrument and systematic noise, reduced flux from 
sky absorption and effects of other varying observing 
conditions that have been omitted. To account in an 
approximate way for these omitted quantities, we re- 
duce the coefficient of the S/N from photon statistics 
and sky noise until S/N matches the empirical data of 
Albrow et al, (1998). For D = 91 cm and t = 300 sec, 
Albrow et al. obtain 1%, 2%, 7% photometry for 
/ = 15, 17, 19 respectively. For / = 19, substitu- 
tion of these values of D and t into Eq. (||) yields 
S/N = 57 instead of a little over 14 for the 7% pho- 
tometry found by Albrow et al. The empirical datum 
is matched if we multiply Eq. (||) by the indicated fac- 
tor of (1/4). Table 1 compares the S/N determined 
by Eq. (|]) with the empirical values of Albrow et al. 
for the remaining two magnitudes and for empirical 
values for a 76 cm telescope at Lick Observatory (W. 



Li and A. Filippenko, Private communication, 2000). 
The Albrow et al. data are nearly matched for all the 
magnitudes, whereas the Lick data are reproduced for 
the brighter magnitudes, but less well for the magni- 
tudes near the limiting value for the telescope. Eq. 
(^J) as modified is thus a good approximation at least 
to the Albrow et al. data, and reproduces much of 
the Lick data. These comparisons with empirical data 
encourage the use of Eq. (^) for our purposes here. 
In determining S/N to be used in finding the prob- 
ability of detecting a planetary companion of a lens 
during an event, we shall assume D = 200 cm and 
t = 60 sec, except we shorten the integration time for 
bright sources to keep S/N < 200. Systematic effects 
would probably soon dominate photon statistics for 
the brighter sources, so S/N would not be improved 
by longer integrations. 



I mag 


Albrow et al. 


Lick 




D=91 cm 
t=300 sec 
s=1.4 as 


D=76 cm 
t=300 sec 
s=2.24 as 


15 


100, 122 


> 80, 100 


17 


50, 46 


35, 35 


18.5 




6, 13 


19 


14, 14 


2-3, 9 



Table 1: Comparison of empirical determinations of 
S/N with those determined by Eq: (^) including the 
factor 1/4. The calculated values of S/N are in bold- 
face type next to the empirical values. 

Subjective impressions of examples of precision 
photometry in the literature, (e.g. Charbonneau et 
al. 2000), suggest that a 2a perturbation of a single 
lens light curve should be easily discerned with a suf- 
ficient number of data points, where a = 1/(S/N) is 
the rms fractional deviation of points about the best 
fit curve (e.g., a = .01 is 1% photometry). What 
is "sufficient" is also subjective, but we can quantify 
this in terms of confidence. We consider a microlens- 
ing event of duration 30 days (2t E ) perturbed by a 
planet with q = 0.0001, where the expected dura- 
tion of the perturbation is ~ 2^/qt E = 0.3 days (e.g., 
Peale, 1997). As an example, we assume a worst 
case true perturbation of the single lens light curve 
remains at the minimum 2/ (S/N) for the entire 0.3 
days, where S/N is that for a single observation of 
the amplified brightness of the source, which will not 
change significantly during the short time of the per- 
turbation. This example perturbation is unlike a true 
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planetary perturbation, but its simplicity and mini- 
mal nature will illustrate the viability of our detection 
criterion adequately. We want to find an estimate of 
the number of data points necessary to detect such 
a perturbation coupled with a measure of the confi- 
dence in the detection. 

In Fig. [l], we show a segment of a microlensing light 
curve with u m in = 0.2 for a source magnitude I = 20 
with this minimum assumed perturbation that occurs 
about 4.5 days after the time of the peak amplifica- 
tion. The solid bold line is the true amplification of 
the source, and the dotted curves span the true un- 
perturbed light curve by a fraction ±2 /(S/N). S/N 
as determined by Eq. (||) has a relatively poor value 
of about 15 at the region of the curve where the per- 
turbation occurs. The simulated data points are nor- 
mally distributed with a for the fractional deviation 
equal to the local value of 1/(S/N). The sampling 
frequency is once per hour outside the perturbation 
region, but it is increased to once every 5 minutes dur- 
ing the perturbation with the same conditions (inte- 
gration time, telescope, etc.) determining S/N as for 
the remaining part of the light curve. The two short 
lines closest to and on either side of the bold line are 
least squares fits to a straight line for the 86 simulated 
data points and to the same set with opposite signed, 
equally probable deviations. The uncertainty Odet i n 
the mean position of the normalized data points dur- 
ing the perturbation is a /\[M = 1/(V86S/N), which 
is a measure of the uncertainty in the detection. The 
effective S/N for the detection is more than 9 times 
that for each individual point and the detection is es- 
sentially certain. That is, with 99% confidence, the 
true amplification is within 3<7det — 1/(3S/N) of the 
measured mean represented by the midpoint of the 
either of the two short lines segments closest to the 
true light curve during the perturbation. The uncer- 
tainty is small compared with the magnitude of the 
fractional perturbation of 2 /(S/N). 

The other two small line segments separated fur- 
ther from the true light curve correspond to least 
squares fits to only 20 points taken during the per- 
turbation, which is a sampling frequency about 3 
times that for the other parts of the light curve. Here 
Vdet — 1/(V20S/N) and there is 99% confidence that 
the true amplification falls within about 2/(3S/N) of 
the midpoint of cither line. This is still sufficiently 
less than the 2 /(S/N) perturbation that a detection 
would be secure. There is a necessary cushion here, 
because the entire light curve will only be a best fit 
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Fig. 1. — Segment of a microlensing light curve for 
/ = 20, u m in — 0.2, and q = 0.0001, where the du- 
ration of the planetary perturbation is 2t E ^fq = 7.2 
hours for an event duration of 30 days. The bold lines 
represent the true amplification and the dotted lines 
are displaced by a fraction of ±2 /(S/N) from the true 
amplification. The simulated points are gaussian dis- 
tributed with a = 1/(S/N) The two short thinner 
lines on either side of and closest to the bold line dur- 
ing the perturbation are least squares fits to the set 
simulated data points and to a set where each devia- 
tion is given the opposite sign. The second set of lines 
further separated from the true amplification during 
the perturbation are least square fits to only 20 data 
points. 

model of all of the points, and the perturbation will 
be suppressed somewhat during the fit. If we assume 
a planetary perturbation of the single lens light curve 
that exceeds 2 /(S/N) can be detected with 20 consec- 
utive data points during the perturbation, confidence 
in the detection would exceed 99% even for the min- 
imum perturbation. A detection criterion can then 
be defined as a fractional perturbation of the single 
lens light curve that exceeds 2 / (S/N) for at least 20 
consecutive data points. Sufficient duration of a per- 
turbation to satisfy this criterion will obtain for al- 
most all of the events considered in determining the 
probabilities for reasonable sampling rates with the 
applicable values of q, and duration will not appear 
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explicitly in the calculations to follow. In this sense, 
duration is really not part of the detection criterion 
used in the calculations, but the constraint of "20 con- 
secutive data points" is understood to be satisfied for 
nearly all detectable events. Note that for the min- 
imum fractional perturbation of 2/ (S/N), some part 
of the perturbation will probably have passed before 
a higher sampling frequency and perhaps longer inte- 
gration times were effected to more closely map the 
perturbed part of the light curve. Hence, a detec- 
tion does not necessarily imply sufficient definition of 
the light curve to constrain q adequately. In addi- 
tion, the poorer definition of the slope of the best fit 
curves when only 20 data points are available illus- 
trates the necessity of very frequent sampling if one 
hopes to characterize the planet with any certainty. 
It is certainly clear from Fig. |l| that an S/N consid- 
erably better than 15 is desirable for characterizing 
the planet, but the smaller values do not preclude de- 
tecting a 2 /(S/N) perturbation with 20 data points. 

The S/N used to establish the criterion for de- 
tecting a planetary companion during a microlensing 
event depends of course on the apparent brightness 
of the source star. The distribution of the apparent 
brightness of the source stars at a particular distance 
are determined from the following LF. 



1.470 x icr 9 e - 930/ * 10.5 < I a < 17 
2.575 x io- 8 e a762/8 17 < L < 18.5 



1.316 x 10~ 4 e 



4^0.300/8 



18.5 < I s < 26, 
(A5) 



which represents the fraction of the sources per unit 
/ magnitude for bulge stars at a distance of 8 kpc. 
Is is the I magnitude at that distance. This /-band 
LF is constructed from the J-band LF of Zoccali et 
ol. (1999) in Appendix A, where the /-band is cho- 
sen as the most likely bandpass that will be used in 
any dedicated planet search from the ground. In Eqs. 
(A5), 10.5 < Is < 17 corresponds to giant stars and 
18.5 < Is < 26 to main sequence stars below the 
turnoff with the second of Eqs. (A5) bridging the 



gap between the two distributions. The same lumi- 
nosity function will be assumed to apply to the disk 
as well as the bulge, although there will be fewer gi- 
ants among the disk stars. Most of the sources will 
be in the bulge in any case, and changes induced in 
the probability calculations by a more refined LF are 
negligibly small. 

For particular LOS, £, q, a and MF of the lenses, 



the probability of detection is calculated for a few tens 
of values of the impact parameter (0 < u m i n < 1) 
and the probabilities averaged over the impact pa- 
rameter distribution. This is done for approximately 
10 different apparent magnitudes in the range / m ; n < 
/ < 21. The maximum value of / = 21 is chosen 
as the dimmest star that can be usefully monitored 
in a ground based search for planetary perturbations, 
whereas / m i n is determined by the brightest star in 
the LF translated to distance £ and increased by the 
extinction to that distance. The extinction is han- 
dled by assuming that the measured extinction to 8 
kpc along the particular LOS is uniformly distributed 
from the observer to 8 kpc. The transformation of the 
LF of Eq. ( |A5| ) to one for the apparent magnitude at 
distance £ ^ 1 is effected by substituting 



h = / apP - A I8 ( - 5 log C 



(5) 



into the LF, where A I8 is the extinction to 8 kpc. 
Details are given in Appendix A. 

3. Probability of detection 

The basic idea behind the detection probability is 
due to Gould and Loeb (1992), although we use an 
S/N criterion for detection instead of a fixed percent- 
age perturbation of the microlensing light curve. Fig. 
^ shows loci of equal amplification for the positions 
of a planet projected onto the lens plane with planet- 
star mass ratio q = 0.001. The amplification of the 
source is A = [1 ± 2/(S/N)]Ao if the planet is on a 
contour, where Aq is the unperturbed amplification of 
the source if the planet were not there. A fractional 
perturbation of 2 /(S/N) is assumed detectable with 
(S/N) being that appropriate to the instantaneous 
brightness of the amplified source. Inside the loci the 
perturbation will exceed 2 /(S/N) in magnitude. The 
method of estimating S/N for a given source magni- 
tude is shown in Section above. 

As the source passes by the lens from left to right 
in Fig. ||, the contours sweep across the lens plane — 
waxing and waning as the source gets closer to or 
further from the lens and as S/N varies with the in- 
stantaneous value of Aq. The large curves trace the 
extremes of the equal amplification loci during the 
traverse of the source across the Einstein ring, and 
a planet would be detected at some time during the 
event if its projected position were anywhere within 
either pair of the large curves. We adopt the proba- 
bility of detecting the planet during the event as the 
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Fig. 2.— Sweep of small contours of ±2/{S/N) 
fractional perturbation of light curve across the lens 
plane. The large curves are traces of the extremes 
of the contours during the event, and they define the 
area within which a planet would be detected if it 
happened to be there. The source magnitude is 1=19. 
The diagram is centered on the lens with the position 
of the source corresponding to the contours displaced 
a small distance into the second quadrant. The dot- 
ted circle is the Einstein ring. At the source loca- 
tion shown, 2/{S/N) = 0.0366, whereas at closest 
approach 2/(S/N) — 0.0283. The source moves hori- 
zontally across the diagram. 



probability that its projected position is within the 
area contained by the large curves. This area expands 
and the probability of detection is correspondingly 
higher as the source brightness is increased thereby 
increasing S/N and as the impact parameter of the 
source-lens encounter is reduced. The contours shrink 
or expand approximately as y/q, and the area be- 
tween the trace curves and probability of detection 
are changed approximately by the same factor. 

Fig. H shows how the area is increased with source 
brightness. The effect on the probability of this 
change in area is shown in Fig. ||, where the detec- 
tion probability is increased for brighter stars for all 
Umin- Fig. H also illustrates the increased probability 
of detection for smaller u m i n . The impact parame- 
ter of 0.07 in Fig. | is reduced to 0.005 in Fig. §, 



which shows the expanded loci of minimum pertur- 
bation for detection and the expanded curves trac- 
ing the extremes of these loci. One notices first in 
Fig. H that the loci extend beyond the traces of their 
extremes along their symmetry axes (the source-lens 
line), so that the probabilities are slightly underes- 
timated by our adopted scheme. That is, a planet 
would still be detected if it were slightly outside the 
curves defining the traces of the extremes of the loci 
contours. However, this effect is most pronounced 
for the smallest impact parameters where it makes at 
most a few percent difference in the calculated prob- 
ability, and we average these probabilities over the 
distribution of impact parameters, which further re- 
duces the relative contribution of this extra probabil- 
ity for small impact parameters. Notice next that the 
contours within which the planet would be detected 
completely encompass the region around the Einstein 
ring radius for small impact parameters. This illus- 
trates the result of Griest and Safizadeh (1998) who 
found unit probability of planet detection for high 
amplification events if the projected position of the 
planet were anywhere within a region spanning the 
Einstein ring radius. 
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Fig. 3. — Traces of the extremes of the contours of 
±2 / (S/N) fractional perturbation of the light curve 
for q = 0.001, Umin = 0.07 and two stellar I- 
magnitudes. The area inside the curves and the prob- 
ability of detection is higher for brighter stars. 
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The amplification of a source by a binary star- 
planet lens is determined in Appendix B with the 
analysis in the complex plane by Witt (1990). The 
procedures for determining the contours of equal am- 
plification and for determining the extremes of these 
contours — used in Figs. |[ ||| and || — are also outlined 
in this appendix. 

4. Probability distribution of the planet in 
the lens plane. 

The time average separation of a planet from its 
primary is a(l + e 2 /2), where a and e are the semi- 
major axis and eccentricity of the orbit. So we can 
account for the random orbit orientation and random 
phase by assuming that the instantaneous position 
of a planet is uniformly distributed over the surface 
of a sphere of this radius. However, we drop the e 2 
and assume the sphere has radius a, where it is un- 
derstood that this "semimajor axis" is augmented by 
some assumed mean orbital eccentricity to represent 
an averaged separation. The cumulative probability 
that a planet with this distribution will have a posi- 
tion projected onto the lens plane between r < a and 
a is (Gould and Loeb, 1992) 

F(r) = f^, 
or between x p = r/R E and a/R E , 



F{x r ,a,(,z,M) = ^l-fii-E 

where K = 4:GM e D s /(c 2 a 2 @ ) = 65.1279 with G = the 
gravitational constant, M Q = the solar mass, a e = 1 
AU, M is the lens mass in units of M Q , and a is now 
measured in AU in the last line. 

We wish to determine a planet detection probabil- 
ity during a microlensing event in terms of a criterion 
based on a representative S/N of the signal. The lat- 
ter depends on the apparent magnitude of the source, 
which in turn depends on £ for the given luminosity 
function. We therefore average the cumulative prob- 
ability in Eq. (0) over the distribution of the lenses 
from the observer to a particular source distance £ 
and simultaneously over the MF of the lenses, since 
the lens mass M is embedded in the definition of R E 
(Eq. (0)). The probability of detection of a planet 
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Fig. 4. — Detection probability for q — 0.001, £ = 
1.0, a — 3 AU as a function of the impact parameter 
for two different source brightnesses. The probability 
increases as u m in is reduced, and it is increased for 
all impact parameters for the brighter star. 

of given q (averaged over the lens density distribution 
and MF about a star lensing a source of given / mag- 
nitude at given distance £ can then be averaged over 
the apparent brightness distribution from the LF as 
adjusted for that particular distance. After the prob- 
abilities are averaged over the MF, lens spatial density 
distribution and LF, they will then be averaged over 
the distribution of visible sources along the particular 
line of sight — all this as a function of the semimajor 
axis of the planetary orbit in AU. But first we must 
average F(x p , a, £, z, M) over the lens distribution be- 
tween observer and source and over the MF. 

The details of this averaging process are given in 
Appendix C. In the average represented by Eq. (pl|), 
a and £ are fixed, and the average is performed as 
a function of x p . The Holtzman et al. (1998) MF, 
which we represent as being oc M^ 1 for 0.08 < M < 
O.7M and oc M~ 2 - 2 for 0.7 < M < 2M , is used 
in the averaging. The lens distribution is determined 
by the triaxial Galactic bulge model of Zhao (1996) 
combined with the double exponential disk model of 
Bahcall and Soneira (1980). 

Fig. H shows the averaged cumulative probability 
F(x p ,a, £) of the projected separation of the planet 
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Fig. 5. — The contours defining a 2/{S/N) pertur- 
bation of the light curve grow and the area inside 
large curves that trace the contour extremes along 
the source-lens line is increased as is the probabil- 
ity for detection for smaller impact parameters u m in- 
The area enclosed by the contours that is beyond the 
traces of the extremes means we underestimate the 
true probability of detection by a few percent for the 
smallest impact parameters. Here the fractional per- 
turbation is fixed at 0.1. 



from its star for several values of the orbit semima- 
jor axis a for the particular choice of a LOS toward 
Baade's window (£,b) = (l°,-4°) (t, b are Galactic 
longitude and latitude respectively) and for a source 
distance of 8 kpc (C = 1). The Holtzman et al. (1998) 
MF is adopted for both disk and bulge stars. Of 
course the cumulative probability is unity for x p = 0, 
since it is certain that the projected position of the 
planet will be somewhere between and a. The only 
contributions to F(x p ,a, £) come from those lenses 
along the LOS where R E < a/x p (Appendix C). For 
small a, the number of lenses whose R E values qualify 
plummets rapidly as x p increases. Increasingly small 
M with z either very small or very near £ are re- 
quired until there are no lenses at all along the LOS 
that satisfy the condition. F(x p , a, C) in Fig. ^ corre- 
spondingly decreases most rapidly toward zero when 
x p increases for small a and less rapidly for larger a. 
Useful in understanding Fig. [| Fig. [?] shows the 



Fig. 6. — Cumulative probability F(x p ,a, C) that the 
projected position of the planet lies between x p and 
a I R E for several values of semimajor axes a for a fixed 
source distance £ = 1. The probability is averaged 
over the line of sight distribution of lenses according 
to the Zhao bulge model and the Bahcall and Soneira 
disk model, and over the Holtzman et al. mass func- 
tion. 



Einstein ring radius averaged along the LOS to dis- 
tance C according to 



(R E ) = (8.07 AU)x 



M„ 



M mi „ 



C .lM Z (C- Z ) z2 dn L (M Z ) dzdM 



dM 



M„ 



c 2 dn L (M,z) 



(7) 



dM 



dzdM 



where dn(M, z)/dM (given in Appendix C) is the 
number density of lenses per unit mass interval at 
position z < C, along the LOS. The integrals are evalu- 
ated with a Monte Carlo technique (Press et al. 1986). 
(R E ) is shown for two MFs in Fig. [7], but concentrate 
on the one for the Holtzman et al. mass function for 
now. The relatively small value of (R E ) for £ = 1 
results from the facts that the MF is dominated by 
M stars with small R E and that R E decreases as the 
lens approaches the source at C = 1 where the num- 
ber density is highest. The curve for the modified 
Holtzman et al. MF (oc M~ 2 ' 2 in Fig. 0) will be 
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used later to illustrate the effect of changing the MF 
on the detection probability. The probability density 



Holtzrnan MF (M~ l for M Stars) 
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Fig. 7. — The Einstein ring radius averaged over the 
lens density distribution along the line of sight toward 
Baade's window (£,b) = (1°,— 4°) and over the indi- 
cated mass functions for the Zhao (1996) bulge and 
Bahcall-Soneira (1980) disk models. 

of the planet separation is the negative derivative of 
these curves so the planet has the highest probabil- 
ity of being found at a value of x p where the slopes 
are steepest. The curves approach zero slope more 
rapidly for the smaller semimajor axes, since it is less 
likely to find a close planet much outside the Einstein 
ring radius for the given MF. Similarly the planet is 
more likely to have a large value of x p if the semima- 
jor axis is larger, so the value of x p where the slope 
is steepest for a given a increases with a. The mag- 
nitudes of the slopes decrease as a increases since the 
probability density is spread over a greater area and 
must be thereby smaller. The probability density of 
the projected position is sharply peaked near r = a 
(Gould and Loeb, 1992). This is understood intu- 
itively since the planet spends a higher fraction of 
the time near the extremes of its projected orbit as it 
travels toward or away from the observer. Or from the 
uniform distribution over the surface of the sphere, a 
large area of the sphere is covered by small changes in 
the LOS near the extreme. The steepest slope for, say, 
a = 2 AU in Fig. § occurs near x p = 1 , which reflects 



the fact that the average R E over the LOS to £ = 1 is 
about 1.7 AU. (That r/R E and r prefers to 

be near a because of the peak at a in the probability 
density.) The steepest slope for a — 10 AU occurs 
near x p — 5 or 6, which is again the planet preferring 
to be at an extreme separation with (R E ) ~ 1.7 AU. 

Fig. U shows the same cumulative probability dis- 
tribution as Fig. |^, but now a is fixed at 3 AU and 
the distance to the source is varied. Again contribu- 
tions to F(x p , a, C) come only from those lenses with 
R E < a/x p . From Fig. we see that (R E ) is large for 
large £, so as x p increases, the number of lenses along 
the LOS contributing to F(x p ,a,C t ) decreases most 
rapidly for large £. That is, if (R E ) is large, the num- 
ber of lenses with sufficiently small R E to contribute 
to F(x p , a, £) plummets most rapidly as x p increases 
for large £, and the number decreases less rapidly for 
smaller £ when x p increases. These trends are evident 
in Fig. ||. 




5 10 15 
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Fig. 8. — Same as Fig. |^ except the cumulative prob- 
ability is given for several values of source distance £ 
for fixed semimajor axis a = 3 AU. 

Since the probability density of the projected planet 
position peaks near r = a, the planet is most likely 
to be separated from the lens by ~ 3 AU for all the 
curves in Fig. ^. The steepest part of the curves in- 
dicating the highest probability density should then 
occur near where x p = 3AU/(i? B ). From Fig. 0, 
(R E ) w 0.8, 2.05 AU at C = 0.2, 1.2 respectively, 
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and we see that the steepest slopes indeed occur near 
x p = 3 w 3AU/0.8AU for C = 0.2 and near x p = 1.5 » 
3AU/2.05AU for C = 1.2. Another way to look at the 
behavior of F(x p , a, C) in Fig. || is to notice that for 
small C, is small so there is a greater probabil- 

ity of finding the planet at larger values of x p . The 
consistencies in the positions of the curves as a func- 
tion of the parameters in both Figs. || and || give one 
some confidence that the cumulative probabilities are 
correctly averaged. 

The following describes how the integral of the 
probability density over the area swept in the lens 
plane out by the contours of 2/(S/N) fractional per- 
turbation of the light curves is carried out using 
F(x p , a, £). This is the integral of dF(x p , a, £)/ 1 dx p in 
Eq. ([|). Assume for the present that the normalized 
distances from the lens to the inner and outer traces of 
the amplification contour extremes in upper half plane 
of Fig. |^ are n/R E and ri/R E respectively. Since 
F(x p ,a, £) is the averaged cumulative probability of 
finding the projected position of the planet in the an- 
nular region between x p and a/R E centered on the 
lens, the probability of finding the planet in a narrow 
wedge of angular width A9 and lying between ri/R E 
and r 2 /R E is just [F(x p — r\j 'R E ,a,C) — F(x p = 
r<2,jR E , a, £)] A9/2ir. The probability of finding the 
planet in the area inside the two trace curves is then 
the sum of these probabilities for all the positions of 
the source as it is stepped through the encounter. A 
similar procedure determines the probability that the 
projected position of the planet lies in the area en- 
closed by the trace curves in the lower half plane of 
Fig ||. The sum of these two probabilities is the prob- 
ability of detecting the planet during this particular 
encounter. This probability is of course dependent 
on the estimate of the S/N (Section ||) we use to de- 
fine the perturbation contours of dctectability in, for 
example, Figs, || and |[ 

5. Procedure 

For each impact parameter, the source is located 
by its angular position 9 measured from the point of 
closest approach to the lens, and the source is stepped 
along its trajectory in equal increments A9. The first 
point is the first encounter of the source with the Ein- 
stein ring of the lens where the unperturbed ampli- 
fication Aq = 1.34, which is here assumed to be the 
minimum amplification for an alert to the event. The 
brightness of the source is increased by the instan- 



taneous amplification to determine the S/N at that 
position. The extremes in the contours correspond- 
ing to the perturbation of 2 /(S/N) are calculated for 
that S/N as outlined above and in Appendix B, and 
the probability that the planet lies in the wedges of 
angular width A9 along the source-lens line and be- 
tween the extreme curves allowing detection is calcu- 
lated for a particular semimajor axis. These "wedge" 
probabilities are summed for the encounter out to a 
source position along its trajectory of 2R E beyond its 
point of closest approach. The detection probability 
is determined for each of a distribution of approxi- 
mately thirty impact parameters for the given source 
distance and magnitude, and averaged over the im- 
pact parameters. 

The calculation is repeated for each of 11 semi- 
major axes (0.6 < a < 10 AU), and one obtains a 
series of functions of a, labeled by the I magnitude 
of the source at the distance (. These are probabili- 
ties P(a, £, I) of detecting the planet during the event 
averaged over the MF of the lenses, over the distri- 
bution of the lenses along the particular LOS from 
the observer to the source at £, and over the distri- 
bution of impact parameters. Figs. ^ and |l^ show 
examples of these curves for £ = 0.2 and 1 respec- 
tively for a selection of source / magnitudes. In these 
figures, q = 0.001, the LOS is toward Baade's win- 
dow (£,b) = (l°,-4°) and the extinction A I8 = 0.76 
mag (Holtzman et al. 1998). At C = 0.2, the bright- 
est star in the LF has 1=7.16, and the curves for 
7.16 < I < 14.5 are superposed, because for these 
bright stars, S/N exceeds 200 for a 60 sec integra- 
tion time, and we have assumed the integration is 
reduced so that S/N does not exceed 200 in such cir- 
cumstances. This means that the detection proba- 
bility is observationally controlled to be the same for 
the given LOS and parameter values for all sources 
with / < 14.5. Otherwise, the detection probabil- 
ity increases with source brightness in these figures 
because the higher S/N increases the area contained 
within the curves tracing the extremes in the ampli- 
fication contours as shown in Fig. |^. 

We pointed out above that the probability distri- 
bution of the projected position of the planet on the 
lens plane is strongly peaked near r = a, and because 
the perturbation of the light curve decreases mono- 
tonically as the planet is moved radially away from 
the Einstein ring radius in either direction (provided 
it is outside the closed critical curves), the probabil- 
ity of detection is maximized when the planet is near 
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Fig. 9. — Probability P(a, C, /) (thin lines) of detect- 
ing a planet when the source at £ = 0.2 for a selection 
of source magnitudes plus the average of the proba- 
bilities over the LF P(a,() (heavy lines). I — 7.16 
is the brightest star in the LF for A IS = 0.76, and 
a single curve applies for / = 7.16 to about 14.5 be- 
cause of our constraint that S/N not exceed 200. The 
probability has been averaged over the LOS distribu- 
tion of lenses, over the Holtzman et al. mass function, 
and over the distribution of impact parameters u m i n , 
where the Zhao bulge model and the Bahcall-Soneira 
disk model are assumed for the stellar mass distribu- 
tion in the Galaxy. 

the Einstein ring of the lens. One then expects the 
probability of detection to be maximized when the 
semimajor axis of the planet is near the average value 
of the Einstein ring radius between the observer and 
the source location. From Fig. Q, at £ = 0.2, 1.0, 
(R E ) « 0.85, 1.7 AU for the Holtzman et al. (1998) 
MF, where relatively small values result from the 
dominance of M stars in the MF. Correspondingly, 
the peak in the detection probabilities occur near (1 
AU, 2 AU) for C = (0.2, 1.0) 

The probabilities are shown for two distributions 
of impact parameters in Figs. ^ and [l^. The random 
nature of the events means that the impact param- 
eters will be uniformly distributed, where we choose 
< u m in < 1 with u = 1 corresponding to an am- 
plification of 1.34 being the assumed threshold for an 



Fig. 10. — Same as Fig. || except the source is at 
C = 1.0 and the brightest star in the LF is / = 11.26. 



alert. However, events with smaller impact param- 
eters are more likely to be observed because of the 
higher magnifications, so the distribution of impact 
parameters of recorded events will be biased toward 
the smaller values. Accordingly, the second distri- 
bution of u min is that of Fig. 9 from Alcock et al. 
(2000), which has a higher proportion of small impact 
parameters to reflect the observational efficiencies. In 
fact, events with u m in > 0.9 contribute negligibly to 
the average of the probabilities over the distribution 
of u m i n , so we choose 30 values 0.01 < u m in < 0.91 
for the uniform distribution. The distribution deter- 
mined by observational efficiencies naturally reduces 
the number of larger impact parameters in the range 
< u min < 1. The higher probabilities for this latter 
distribution of u m i n reflect the higher probabilities for 
closer impact parameters implied in Figs. || and ^|. 

Also shown in Figs. || and [h] are the averages 
over the distribution of / magnitudes of the sources 
at each distance. These are effectively averages over 
the LF expressed in / magnitudes appropriate to each 
distance but truncated to the visible sources (/ < 21). 
The averaged probability is 

,21 

P(a,Q = J i P(a,<:,I)h(I)dI, (8) 
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where I min = 10.5 + A J8 C + 51ogC and f L (I) is f L {h) 
in Eq. JA5| ) with Eq. (||) replacing I s and the result 
renormalized so that f L (I)dI represents the fraction 
of visible sources in range dl about I. Examples of 
P(a, C, I) are the labeled curves in Figs. || and and 
P(a, C) is labeled with (I). Since the curves are for 
discrete values of /, the integral in Eq. (||) is approx- 
imated by J2i P i a >C,Ii)fL(Ii)AIi, where h are the 
discrete values of / at which P is evaluated and AJ^ 
is an increment centered on 7, and stretching halfway 
to adjacent values of Ii on cither side. Fig. [ll] demon- 
strates this procedure, where the shaded area under 
the histogram is the average over the LF for ( = 1 
and a = 2.2 AU. Fig. |l2| shows the detection prob- 
ability averaged over the luminosity function of the 
sources for the uniform distribution of impact param- 
eters, over the Holtzman et al. MF and over the LOS 
distribution of lenses for a series of values of The 
semimajor axis where the probability is maximal fol- 
lows the mean value of R E as described above. The 
derivation of the overall averaged detection probabil- 
ity will be completed by the average over the distri- 
bution of visible sources. 

6. Average of P(a, () over( 

The final average of the detection probability is 
carried out with the following expression: 

f P(a,0n s (0t 2 ttF2i,C,A IS )d( 
P(a) = ^-o (9) 

/ n s (oc 2 e(^2i,c,^R 

Ja.i 

where n s (C) is the number density of sources at dis- 
tance C and where 

Z(F ,CA IS ) = C 1 {AnDlT 21 C 2 e- omiA '»Y+C 2 , (Cll) 

derived in Appendix A, is the fraction of the sources 
with observed flux densities T greater than a min- 
imum T21 assumed detectable. T21 — 6.21 x 10~ 15 
ergs cm _2 sec _1 is the minimum flux density in /-band 
considered observable — here from a star with I = 21. 
C\lP la + C2 is the fraction of the sources with lumi- 
nosity L 1 > L I0 , which follows from the fraction of 
sources with I < Iq found from the LF in Appendix 
A. The constants and index are different in the LF for 
each of three ranges of Ig (Eq. (|A5|)) with 1.16% of 
the stars having Jg < 17 and 4.18% of the stars hav- 
ing ig < 18.5, with I = 17 and 18.5 being the values 
where the LF index changes. Hence, Ci, C2 and j3 
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Fig. 11. — An example of the average of the detection 
probabilities over the LF for ( — 1 and a = 2.2 AU. 
The probabilities have already been averaged over the 
the distribution of lenses along the LOS, over the 
Holtzman et al. mass function and over a uniform dis- 
tribution of impact parameters. The shaded area un- 
der the histogram is the averaged probability, which 
is 0.113 in this case. 

must change when £ = 0.0116 and 0.0418. Also £, = 1 
for values of £ so small that all of the sources are vis- 
ible and £ = for values of £ so large that even the 
brightest source in the LF has I > 21. The parame- 
ters in £ and examples of values of C where changes 
occur for a few values of A I8 are given in Appendix 
A. Fig. |l3| shows £ as a function of £ for three values 
of A I8 . 

We choose source distances for which the detection 
probability is to be determined to range from ( = 0.1 
to 1.2 in increments of 0.1. It is very unlikely for any 
source to be closer to us than £ = 0.1 (0.8 kpc), and 
beyond C — 1-2, extinction and blending will proba- 
bly make the number of sources that are usably visible 
negligibly small. The source number density along the 
LOS is obtained from the Zhao (1996) bulge model 
and the Bahcall-Soneira disk for the Holtzman et al. 
(1998) mass function in Appendix B. The overall co- 
efficients in the bulge and disk distributions are de- 
termined by the constraints on the total bulge mass 
of 2.2 x 1O 1O M (Zhao, 1996) and a local disk mass 
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Fig. 12. — Average of the detection probability over 
the LOS distribution of lenses, over the Holtzman et 
al. MF, over a uniform distribution of impact pa- 
rameters Umin, and over the LF for £ incremented in 
steps of 0.1 (except £ = 0.4 is not included). The 
peaks occur where a « (R E )- 

density of 0.05M o pc -3 (Bahcall and Soneira, 1980) 
respectively. The functions P(a,() (E.g., Fig. |l^) 
are made continuous functions of ( by interpolating 
between the discrete points where they are evaluated. 
This makes the integrands for a given LOS in Eq. (^) 
functions of £ alone and evaluation of the integrals is 
straightforward. 

Fig. |lj shows the averaged probability P(a) for 
a uniform distribution of impact parameters and for 
that distribution appropriate to the efficiency of de- 
tection by the MACHO group (Alcock et al, 2000). 
The line of sight is toward Baade's window (I, b) = 
(l°,-4°), where A IS = 0.76 (Holtzman et al. 1998). 
The higher weight given to the small impact parame- 
ters in the latter distribution nearly doubles the prob- 
ability of detection obtained from a uniform distribu- 
tion of u m in- The large probabilities at small planet 
separations that resulted for small ( (Fig. |l2|) are seen 
to be suppressed by the lack of sources at small (. The 
semimajor axis for the peak probability of detection 
is near 1.5 AU, which is close to the mean R E over 
the entire LOS for the Holtzman, et al. mass function 
and Zhao plus Bahcall-Soneira galactic model. 



Fig. 13. — Fraction of sources visible for the adopted 
LF as a function of source distance. F21 is the flux 
density from a star with /-magnitude = 21 that is 
assumed to be the minimum observable, and A I8 is 
the /-band extinction to 8 kpc. The different line 
types correspond to the different forms of the LF for 
different brightness ranges. 

7. Parameter dependence of the detection 
probability 

Fig. [l5| shows the reduction in the probability 
by increasing the extinction to A I8 = 3, by assum- 
ing a 3<7 instead of a 2a detection both along the 
LOS to Baade's window and by changing the LOS to 
(£, b) = (10°, -4°) with A I8 restored to 0.76 mag. All 
assume the distribution of u m i n appropriate to the 
efficiency of observing the event, and this curve for 
Baade's window and A IS = 0.76 is repeated from Fig. 
|l4| for comparison. The 3/(S/N) detection criterion 
reduces the probability of detection for all semimajor 
axes over that for a 2/(S/N) detection criterion as ex- 
pected, and the mechanism is a reduction in the area 
in the lens plane inside of which a planet's projected 
position would lead to a detection. Changing the line 
of sight to {£, b) = (10°, -4°) (also with the i/{S/N) 
detection criterion) from (£, b) = (1°, —4°) with equal 
extinctions increases the probability of detection for 
a > its value corresponding to the peak probabil- 
ity and decreases the probability for a < its value 
near the peak. By swinging the LOS away from the 
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Fig. 14. — Averaged probability of detecting a planet 
during a microlensing event toward Baade's window 
for a uniform distribution of impact parameters and 
a distribution determined by the observational effi- 
ciency of the MACHO group (Alcock et al. 2000). 



bulge, we have effectively increased the mean Einstein 
ring radii of the lenses. There are considerably fewer 
sources and fewer lenses at the most distant parts of 
the LOS at the larger Galactic longitudes, so there is 
less weight given to those configurations of small Ein- 
stein ring radii in the averages over both the lens and 
source distributions. That means that planets with 
small semimajor axes (near R E where the probability 
is large) will have a smaller probability of detection, 
whereas those with larger semimajor axes will have 
their averaged probabilities increased. One can see 
from Fig. |l5| that this is not a large effect. The in- 
crease in A IS toward Baade's window has the effect of 
increasing the apparent magnitudes of all the sources, 
thereby reducing the S/N and the probability of de- 
tecting a planet. 

Fig. [l6| emphasizes the effect of changing the line 
of sight relative to the center of the Galaxy. Here 
a series of LOS with different Galactic latitudes and 
the same Galactic longitude (£ = 1°) again show the 
effect of decreasing or increasing the proportion of 
the lenses for large distances from the observer. The 
probability of detection is raised for larger semima- 



Fig. 15. — Averaged probability of detecting a planet 
during a microlensing event showing effect of shifting 
the LOS, of increasing the extinction, and of using 
a 3/ (S/N) instead of a 2/(S/N) detection criterion. 
The probability is averaged over the source-lens im- 
pact parameter distribution determined by the obser- 
vational efficiency, over the lens distribution along the 
LOS, over the MF, over the LF, and over the visible 
source distribution along the LOS. 

jor axes for the LOS directed further from the high 
spatial density region near the Galactic center, and it 
is lowered for the LOS closer to the Galactic center 
for reasons discussed above. Rather than show the 
complete averaged probability for these examples, we 
notice that the probability averaged only over the dis- 
tribution of impact parameters, over the distribution 
of lenses along the line of sight to ( = 1.0 and over 
the MF for / = 19.5 mag gives a good approximation 
to the detection probability with additional averages 
over the LF and the distributions of sources leading 
to P(a). The points indicated by stars in Fig. |l6| are 
appropriate to P(a) and indicate how well the partial 
average for the selected values of C and / approxi- 
mates the total average for (I, b) = (1°, —4°). 

The peak probability, in addition to being shifted 
slightly to larger semimajor axes as \b\ is increased, is 
also increased slightly. This is an extension of the in- 
crease in probability for larger semimajor axes as less 
weight is given to the more distant lenses with in- 
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Fig. 16. — Averaged probability of detecting a planet 
during a microlensing event showing the effect of shift- 
ing the line of sight along a line of constant galac- 
tic longitude t = 1°. Here P(a,£,I), with C = 1 
and / — 19.5, is used to approximate P(a), where 
the galactic coordinates have been added in the argu- 
ments to distinguish the LOSs. The stars correspond 
to P(a) = P(a, 1°,— 4°) and indicate the good ap- 
proximation by the partial average for the selected / 
magnitude for(O) = (l°»-4°). 

creasing \b\. Since the averaged probability depends 
in an intricate way on the weighting of the Einstein 
ring radius by the distribution of lenses along the LOS 
relative to the semimajor axis of the planetary orbit, 
there is no reason to expect the peak probability not 
to change as it shifts to larger semimajor axes. There 
is almost no difference in the probability in Fig. [l6| 
for b = -15° and -20°. The LOS for these lati- 
tudes does not intercept much of the bulge distribu- 
tion of stars, so the distribution of lenses along the 
LOS, being mostly disk stars, changes very little be- 
tween b = 15° and 20°. 

The reason for considering the modified Holtzman 
et al. mass function in Fig. ^ is that more M stars 
than are observed seem to be required in the MF to 
account for the relatively large number of short time 
scale events toward the Galactic center. The mod- 
ified Holtzman mass function cx M~ 2 2 for all stars 
0.08 < M < 2.0M o gives a reasonably good ap- 



proximation to the time scale frequency distribution 
obtained by the MACHO group for the 1993 bulge 
season (Peale, 1998), but extending the MF into the 
brown dwarf region does not (Peale, 1999). The larger 
number of small mass stars in this distribution leads 
to the reduction in the averaged Einstein ring radius 
for all £ in Fig. fjj. Its consideration here allows us 
the opportunity to illustrate the effect of changing the 
mass function on the probability of planet detection. 
We expect this reduction in the mean Einstein ring ra- 
dius to shift the peak in the probability distribution 
to smaller semimajor axes, to lower the probability 
for the larger semimajor axes, and to raise it for the 
smallest semimajor axes as discussed above. This ex- 
pected effect is shown in Fig. 17, where we have again 



used the partial average of the detection probability 
as an approximation to P(a). The peak probabil- 
ity has moved from something near 1.8 AU to being 
near 1.4 AU, again showing that the peak occurs near 
the averaged Einstein ring radius (R E ) ~ 1.4 for the 
largest values of £ in Fig. ^ (lower curve). 
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Fig. 17. — Effect on the averaged detection probabil- 
ity of changing the MF. The other information is the 
same as in Fig. |l^. The larger number of M stars in 
the MF with the larger index leads to a shift in the 
peak probability to smaller a because of the smaller 
(R E )- This same shift in (R E to smaller values leads 
to a decreased probability for large a and an increased 
probability for small a. 
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The averaged probability of planet detection varies 
approximately as ^fq. This is shown in Fig. [l8|, where 
decreasing q by a factor of 10 leads to a reduction in 
the probability of about ylO. The actual ratio of the 
probabilities for the two values of q starts at 3.14 for 
a = 0.6, drops to 2.7 near the peak and rises to 5.7 
for a = 10, which shows that the approximate scaling 
is rather crude for some values of a, but it still gives 
a useful extension of the results for q = 0.001 and 
0.0001 shown in Fig. [l8]. This scaling should not be 
carried to mass ratios less than 10~ 4 , since effects of 
the non-zero size of the source become important and 
are not included in this analysis (e.g., Bennett and 
Rhie, 1996). On the other hand, the probabilities 
so obtained may not be that far from the Bennett 
and Rhie values for small planets (Peale, 1997). The 
source of this scaling can be traced to the scaling of 
the size of the contours of equal perturbation in Figs. 
H and H approximately as ^fq (Gould and Loeb, 1992). 
The length of the arc described by the center of the 
contours during the event is independent of q. So 
only the dimension of the total area swept out by the 
contours that is perpendicular to this arc is affected 
by q. As this dimension scales approximately as ^/q, 
the whole area swept out scales approximately as ^fq. 
Since the probability per unit area of the projected 
position of the planet falling into the area usually does 
not vary much over the area, the total probability of 
finding the planet in the area scales roughly as the 
area and hence ~ ^fq. 

We have averaged over the LF in the above prob- 
abilities, but since the magnitude of the source is 
known in any particular event, we could eliminate 
this average and evaluate the detection probability 
for particular / magnitudes. These probabilities are 
shown in Fig. [ll| The probabilities for fixed / mag- 
nitudes averaged over the distribution of lenses along 
the LOS, over the lens MF, over the distribution of 
impact parameters (here determined by observational 
efficiency) , and over the distribution of visible sources 
along the LOS are higher for brighter stars because 
of the increasing S/N with decreasing /-band mag- 
nitude. As in the other examples, the probabilities 
are significant even for the faintest stars that we have 
assumed suitable for followup photometry from the 
ground. 

The peak amplification during an event will yield 
the impact parameter with an accuracy determined 
by how well the blending can be modeled, so the av- 
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Fig. 18. — Effect on the completely averaged detec- 
tion probability of changing the planet-star mass ratio 
q showing that P(a) ~ ^fq. The actual ratio between 
the two curves is 3.14 at a = 0.6, decreases to 2.7 
near the peak and increases to 5.7 at a = 10.0, so the 
scaling is poorest at the larger values of a. 



events. Fig. 20 shows the detection probabilities for 



erage over u Tl 



can also be eliminated for individual 



individual events involving a 19th magnitude source 
with impact parameters u m in — 0.1 and 0.3. The 
probabilities are now averaged only over the degener- 
ate parameters of the distances to the lens and source, 
weighted by the distribution of lenses and sources 
along the LOS, and the MF. The rapid drop in the 
detection probability as the impact parameter is in- 
creased is evident, and it again illustrates how the bias 
toward small impact parameters in the observational 
record increases the completely averaged probability 
of detection in Fig. [l4|. 

Since Fig. ^ is based on idealized observational 
conditions, it is probably of limited use in interpret- 
ing real microlensing data from an individual event. 
However, the techniques developed here can be used 
in the interpretation of the Gaudi and Sackett (2000) 
analysis of null results from real data (Albrow et al. 
2000b) by eliminating the guesses at representative 
source and lens distances. The x 2 analysis of Gaudi 
and Sackett outlined in Section 1 yields a best esti- 
mate of the impact parameter and source / magni- 
tude, which we shall assume as known for the prob- 
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Fig. 19. — Averaged detection probability as a func- 
tion of the I magnitude of the source for the LOS 
toward Baade's window. The average over the LF 
is bypassed here and P{a) — > P(a,I) on the verti- 
cal axis. The impact parameter u m i n is distributed 
according to the efficiency of event detection. The 
Holtzman et al. (1998) mass function is assumed. 

ability e(x p ,q) in Eq. (g). Like Fig this elim- 
inates averages over u m m and over the LF of the 
sources in Eq. (JsJ) . The elimination of the LF av- 
erage means the integration over the distribution of 
sources can be included with the integration over the 
distribution of lenses and over the MF. The point 
where q enters the averaging process is in the inte- 
gration of dF(xp,aX)/dx p in Eq. (||) over the area 
S(I, u l min ) inside of which a planet would be detected, 
since q, I, u m i n determine the extent of the area. To 
determine the probability of detecting a planet with 
a particular mass ratio q at semimajor axis a for the 
real data from a particular event, e(x p , q) is integrated 
over the probability density —dF(x pi a)/dx p over the 
range of x p where the density is non-negligible and 
for which e(x p , q) has been determined from the data. 
Note that £ is missing from the arguments of F(x p , a) 
since we have also integrated over the distribution of 
source distances weighted by the source distribution. 
Hence, 

f 1 ?™. / d \ 

P'ia, 0) = J dx p e(x p , q) ^ _ ^rj 



Fig. 20. — Probability of detecting a planet during 
an individual event where both the impact parameter 
and the / magnitude of the source are known. The 
averages over the impact parameter distribution and 
over the LF are omitted, leaving only averages over 
the degenerate parameters, z, £ and M, where the 
former two are weighted by the distribution of lenses 
and sources along the LOS and the last by the MF. 
P(a) — > P(a, I,u m i n ) on the vertical axis. Note the 
rapid drop in the detection probability as the impact 
parameter is increased. 

,4.2 K fM ma x 

x / / / n' s (0K(z,M) 

JO.l JO JMmin 

x Q[f(x p , a, C, z, M)]F(x p , a, £, z, M)d(, dzdM, 

(10) 

is the probability of detecting a planet with mass ratio 
q at semimajor axis a for the particular event whose 
data were used in the \ 2 analysis to yield e{x p ,q). 
The variables in Eq. ([ll]) are defined just below Eq. 
(|J). The integral is repeated for a sufficient number 
of values of a to define the functional form of P' (a, q) . 
The integrals are to be evaluated with a Monte Carlo 
technique as before. It should be noted that e(x p ,q) 
depends only on the data and is model independent, 
whereas P'(a,q) depends on MF, and the lens and 
source distances weighted by the lens and source dis- 
tributions along the LOS, and finally P(a, q = 0.001) 
in Eq. (||) depends on all of these parameters plus 
the LF of the sources and the distribution of impact 
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parameters. 

8. Summary 

We have derived the probability of detecting a 
planetary companion of a star acting as a gravita- 
tional lens of a more distant star toward the Galactic 
bulge sometime during the microlensing event. The 
probability is determined as a function of the or- 
bital semimajor axis of the planet, and it is based 
on a criterion which assumes that a perturbation to 
the microlensing light curve due to the presence of 
the planetary companion that exceeds 2/(S/N) for at 
least 20 consecutive photometric data points is de- 
tectable. By the observer's choice under the assump- 
tion that many events are being monitored simulta- 
neously and the likely dominance of systematic noise 
for brighter stars, integration times are assumed to 
be terminated when S/N exceeds 200. Given this cri- 
terion, the probability for q = 0.001 (or 0.0001) is 
averaged over the distribution of the projected posi- 
tion of the planet relative to the lens, over the spatial 
distribution of the lenses along the LOS for a model 
of the Galaxy incorporating the Zhao (1996) and the 
Bahcall-Soneira (1980) bulge and disk model respec- 
tively, over the Holtzman et al. (1998) mass function, 
over the distribution of source-lens impact parameters 
inside the Einstein ring radius R E , over the /-band LF 
derived from the J-band LF of Zocalli et al. (2000) 
as adjusted by the source distance and the extinction 
thereto, and over the spatial distribution of sources 
along the LOS, which is identical to the distribution 
of lenses except only a fraction are visible due to dis- 
tance and extinction. 

The averaged probability scales approximately as 
y/q (Fig. |l^), but the scaling should not be extended 
to values of q < 10~ 4 where the effect of the non-zero 
size of the sources becomes important. The detection 
probability has its maximum value of a little over 10% 
if a uniform distribution of impact parameters is as- 
sumed and over 20% if the distribution of impact pa- 
rameters is adjusted for the observing efficiency found 
by Alcock et al. (2000) (Fig. f[§). The maximum oc- 
curs for a planet semimajor axis that is close to the 
average value of R E over the LOS — 1.5 to 2 AU for a 
LOS toward Baade's window {I, b) = (1°, -4°) for the 
Holtzman et al. (1998) MF. The averaged probability 
of detection is decreased with increasing extinction, 
because of the lowered S/N for all sources. Requir- 
ing a detectable perturbation of the light curve to 



be 3/ (S/N) rather than 2/ (S/N) decreases the prob- 
ability by about 20% (Fig. [D|. If the probability 
is not averaged over the LF, the otherwise averaged 
probability increases with source brightness because 
of the higher S/N (Fig. |l|). A MF with a higher 
proportion of M stars shifts the peak in the probabil- 
ity to smaller semimajor axes as expected (Fig. |l7| ). 
It also reduces the probability of detecting a planet 
at large semimajor axes while increasing the proba- 
bility of detecting planets at very small semimajor 
axes. These trends are a consequence of the reduced 
(R E ), where preference for semimajor axes compara- 
ble to R E result from highest probability of detecting 
a planet occurring when its projected position is near 
R E . These same trends in the probability are seen 
when the LOS is moved away or toward regions of 
high source density near the Galactic bulge (Figs. [l6| 
and |l5|) , where the cause is again the change in (R E ) 
because of the different weightings of the source and 
lens distributions in the averages for different LOSs. 

For individual events, good estimates of the im- 
pact parameter and the magnitude of the source are 
usually available. The averages over the impact pa- 
rameter distribution and over the LF are then omit- 
ted in determining the probability of detection for a 
particular event leaving averages only over the degen- 
erate parameters. The drop in detection probability 
with impact parameter for a 19th magnitude source 
is illustrated in (Fig. |2(]). For real individual events, 
the empirically determined, model independent prob- 
ability e(x p ,q) is averaged over probability density 
—dF(x v ,a)/dx v to yield P'(a,q), the probability of 
detecting a planet with given planet-star mass ratio 
q at semimajor axis a for that particular event. 

9. Discussion 

All of the variations of parameters we have inves- 
tigated maintain the relatively high detection prob- 
abilities over a substantial range of semimajor axes, 
while shifting the peak slightly and altering the rela- 
tive probabilities of detecting planets with large and 
small semimajor axes. This relative robustness rein- 
forces the optimism for planet detection through mi- 
crolensing generated by the Gould and Loeb (1992) 
paper. Continuous coverage of the light curve is as- 
sumed, whereas weather and hardware maintenance 
limit the Galactic bulge coverage by even a global 
distribution of 3 or 4 telescopes to a maximum near 
80% at the midpoint of the bulge season, and this 
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falls off on cither side of mid season because a smaller 
fraction of the night is usable (Pealc, 1997). Still, 
for Jupiter and Saturn mass planets, defined in terms 
of the planet-star mass ratio, the perturbation will 
last on the order of one day and would usually be 
detected, although the light curve during the pertur- 
bation would not be completely covered. If breaks in 
the coverage exceed one or two days at a time, the 
probability of planet detection will be reduced by the 
percentage of total event time that is missed. Our 
detection criterion for a perturbation of magnitude 
2/(S/N) persisting for at least 20 consecutive data 
points means that missed coverage during the pertur- 
bation is serious for detection and much more so for 
determining q. 

Necessarily, we have simplified the calculation of 
S/N for unknown observing systems. The S/N from 
photon statistics and sky background is reduced by 
a factor 1/4 to account for systematic and other ne- 
glected noise sources. The factor was chosen to yield 
the empirical S/N values for two real systems with 
their observing conditions. Although this yields rea- 
sonable S/N estimates for our purposes, each sys- 
tem must be separately calibrated, and this constrains 
other parameters such as integration time for desired 
S/N w 100 or more. Increased integration time may 
compromise the desired sampling time of a data point 
every one or two hours if there are too many ongo- 
ing events that must be monitored. Either the overall 
probability of planet detection would be thereby re- 
duced, or the number of events that can be adequately 
followed with high S/N will be limited in a ground 
based search. The averaged detection probabilities 
thus explicitly depend on our assumptions of a 2 m 
telescope, 1.4 arcsec seeing, and a fixed integration 
time of 60 seconds, because of their determination of 
S/N for a given magnitude source. 

We have neglected blended light under the assump- 
tion that image subtraction techniques (Alcock et al. 
2000) can eliminate this contamination for a wide 
range of source magnitudes and lensing magnifica- 
tions. If unmodeled blended light remains, it will 
have the effect of decreasing the fractional perturba- 
tion of the light curve by the planet and reducing the 
probability of detection in addition to introducing an 
uncertainty in the impact parameter. Several param- 
eters, such as integration time, are not going to be 
fixed as we have assumed, but will vary even during 
a single night as conditions change. We have also 
assumed the single lens light curve is known a priori, 



when only a light curve that is the best fit to the data 
will be available. 

The complete averages of the planet detection 
probabilities we have derived, in addition to estab- 
lishing the robustness of their relatively large values, 
are useful in determining the effect of various distribu- 
tions of parameter values. For establishing the mean- 
ing of null results, however, all the information avail- 
able for each particular event should be used. Aver- 
aged values should be used only for those properties of 
the lens-source system that cannot be lifted from the 
degeneracies. The magnitude of the source is going to 
be more or less known, so averages over the LF are not 
appropriate in analyzing individual events. The peak 
amplification will yield the impact parameter with an 
accuracy determined by how well the blending can 
be modeled, so the average over u m i n should also be 
eliminated for individual events. We have shown how 
the probability varies for various magnitude sources, 
but averaged over all other parameters in Fig. [l9|, and 
Fig. ^0| shows how the probability decreases for in- 
dividual events as the impact parameter is increased. 
This decrease can be also be estimated from Fig. |[ 
The degeneracies remaining are the distances to the 
lens and source and the mass of the lens, so for most 
events, averages over the distribution of lenses and 
sources along the LOS and over the MF are still ap- 
propriate. One or more of these averages might be 
removed for individual events if the parallax of the 
lens can be determined. 

There have been several comparisons of the mi- 
crolensing technique for the detection of planets with 
the very successful radial velocity technique (e.g., Al- 
brow et al. 2000b), where the emphasis has been on 
their complementary detection of distant and close 
planets respectively or microlensing's unique sensi- 
tivity to terrestrial mass planets. Microlensing's ad- 
vantage for detecting large, distant planets from the 
ground and compiling meaningful statistics therefrom 
decreases as the time span for the radial velocity pro- 
grams with the highest precision increases. However, 
most of the lenses will be M stars, whereas these stars 
comprise only a small fraction of the radial velocity 
targets. In addition, the M stars will be located pri- 
marily in the Galactic bulge in contrast to the local 
nature of the stars in radial velocity searches. In that 
sense, microlensing searches from the ground will al- 
ways have a unique contribution to make in constrain- 
ing the statistical distribution of planets. 

Its unique sensitivity to Earth mass planets may 
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be less important because of the reduced probabil- 
ity of detecting such planets from the ground. This 
importance would be drastically increased with the 
promise of an enormous increase in the number of 
events and the continuous monitoring of all the events 
with high time resolution, high precision photome- 
try from space (Bennett and Rhie, 2000). Such a 
space based planet detection scheme using microlens- 
ing, if realized, would decrease the motivation for 
ground based searches because of the drastically in- 
creased throughput and higher sensitivity. However, 
no such mission has yet been funded, so it is pre- 
mature to abandon continued planning for extensive 
ground based microlensing searches for planets. The 
robust nature of the large probabilities of detection 
for Jupiter to perhaps Uranus mass planets during 
an event demonstrated here supports the continua- 
tion of this effort. The need for the statistics of the 
occurrence of planets distant from their sources that 
can be obtained relatively rapidly from the ground 
only with a dedicated microlensing search will become 
more urgent as planning for various space-based plan- 
etary searches continues. 
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A. The luminosity function 



= 6.096 x 10~ 5 e - 366 J , 18 < J < 24 (Al) 

represents the fraction of sources per unit J magni- 
tude. Although the extinction is lower in J, the sky 
noise may be too excessive in this band for either the 
initial microlensing survey or the high time resolution 
followup monitoring in search of planetary perturba- 
tions (C. Stubbs, private communication 1999). We 
will therefore convert the LF to /-band, for which 
we need I — J as a function of J. For this purpose, 
we use Tables II and III of Bessell and Brett (1988), 
which give the colors of dwarfs (Class V) and giants 
(Class III) respectively as a function of spectral type 
in terms of magnitude differences. 

V — J and I— J are determined as functions of M v 
with appropriate relations among the columns in the 
tables, and M v is substituted for the spectral types 
separately for dwarfs and giants. M v values are taken 
from § 96 of Allen (1973). Linear fits to the points in 
the V — J, M v and I — J, M v planes yield 

V — J = 0.361 M v - 0.422 for dwarfs, 

= -1.833 M v + 2.667 for giants, (A2) 

I — J = 0.141 M v - 0.159 for dwarfs, 

= -0.556 M v + 0.944 for giants. (A3) 

Since V — J = M v — M 3 we have M v as a function of 
Mj which when substituted into the I — J relations 
with Mj = J — 14.5, where 14.5 is the distance mod- 
ulus for the center of the Galaxy, yields I - J as a 
function of J: 

I — J = 0.221 J - 3.456 for dwarfs, 

= -0.196 J + 3.263 for giants. (A4) 

We solve each of Eqs. (|A|) for J (I). Then the 
LF in /-band is given by f L (I) = f L [J(I)}{dJ/dI). 
However, the range of / for each segment differs 
from the range in J; J = 9, 17, 18, 24 — > / = 
10.60, 16.93, 18.52, 25.85. We have assumed the ex- 
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Zoccali et al. (2000) assemble a luminosity func- 
tion (LF) for stars in the Galactic bulge in J-band 
from their own data for stars later than the main 
sequence turnoff and from Tiede et al. (1995) and 
Frogel & Whitford (1987) for the giants. We approx- 
imate the LF in their Fig. 10 with three straight line 
segments and further normalize the integral over the 
LF to unity such that 

f L (J) = 2.628 x 10" 8 e - 748 J , 9 < J < 17, 
= 1.108 x 10 _1 V- 612 J , 17 < J < 18, 



pressions for the giants apply for 9 < J < 17 and 
expressions for the dwarfs apply for 18 < J < 24. 
round the values of / corresponding to the segn 
boundaries in J slightly and assume the giant i 
tions apply for 10.5 < / < 17 and the dwarf relat 
for 18.5 < I < 26. The giant and dwarf segments 
connected with a straight line in the /-band equ 
lent of Fig. 12 of Zoccali et al. (1999) and the 
renormalized so that 

f L {I) = 1.470 x 10- 9 e a930/s , 10.5 < / 8 < 17 
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2.575 x l(T 8 e a762/8 , 17 < J s < 18.5 
1.316 x 10- 4 e a300/8 , 



18.5 < h < 26(A5) 



represents the fraction of the sources per unit / mag- 
nitude. (The subscript stands for the 8 kpc distance 
to the Galactic center.) The slope of the /-band lu- 
minosity function of Holtzman et al. (1998) for the 
dim stars is in reasonable agreement with the last of 
Eqs. (|A|). 

The luminosity function will be used in two ways 
in the determination of planetary detection prob- 
abilities: First in the S/N criterion for detection 
and second in the visibility of the sources. Both 
will involve a modification of the luminosity func- 
tion to account for extinction as a function of the 
distance to the source. The extinction A IS to 8 kpc 
is about 0.76 magnitudes in Baade's window (Holtz- 
man et al. 1998) but increases to about 18 to- 
ward the Galactic center. A IS will be retained as 
a free parameter to be adjusted for each LOS. We 
assume the extinction is uniformly distributed along 
the LOS such that the fractional change in inten- 
sity dl/l = -CdD os and J = J e- CDos , where 
D os is the distance to the source and C is con- 
stant. Then e - CD ° s = 10-° AA ' = C-0-921A, yieldg 
CD os = OmiAj = CD S C = 0.921A J8 C or the extinc- 
tion = A IS C,. The apparent magnitude of a star 
at arbitrary £ along a LOS where the extinction to 8 
kpc is A I8 is then 



7 = / 8 + yl J8 C + 51ogC. 



(A6) 



For the visibility of the sources it is convenient 
to express the luminosity function in terms of lumi- 
nosity instead of magnitudes, so we proceed to de- 
termine the /-band luminosities of the bulge stars. 
The /-band is centered near 0.80 ^tm, with a flux 
density at the Earth from a magnitude star like 
Vega of 1.20 x 10" 12 Wcm- 2 ^m _1 (McLean, 1997). 
The Gunn / filter spans 0.75 to 0.88 fim (FWH M) 



( http://ftp.noao.edu/kpno/filtcrs/4inchJist.htm] ), so 
the flux density in / band from Vega is 1.56 x 
10 _6 ergs cm _2 sec _1 . Then, 



3.06 x 10 6 10-°' 4/s = 3.06 x io 6 e ~ a921/s (A7) 



is the luminosity of a star at the Galactic center whose 
/ magnitude is 1$. 

The fraction of stars with L I > L I0 is the fraction 
of stars with / 8 < J g0 , which is found by integrating 



Eq. (A5) from the minimum Is value of 10.5 to J g0 . 
The latter fraction is 

F(I 8 < J 80 ) = 1.581 x io- 9 e - 930/so - 2.75 x 10~ 5 

10.5 < / 80 < 17 
= 3.381 x icr 8 e - 762/ so _ 2.590 x 10~ 3 

17 < / 80 < 18.5 
= 4.387 x icr 4 e - 300/ 8o - 0.0711 

18.5 < / 80 < 26 (A8) 



Using Eq. ( JA7| ) to convert /§ to L l7 we find the frac- 
tion of stars with L 1 > L I0 to be 

F{L 1 > L I0 ) = 5.64 x 10" 3 L7 1 01 - 2.75 x 10~ 5 , 
194 > L I0 > 0.487, 
= 7.84 x 10~ 3 L7 ' 827 - 2.59 x 10~ 3 , 
0.487 > L I0 > 0.123, 



5.71 x 10-^L 



-0.326 



- 0.0711, 



0.123 > L m > 1.23 x 10~ 4 , (A9) 

where L t and L I0 are in units of L e . 
Source Visibility 

We mo dify the cumulative luminosity function of 
Eqs. (A9) with extinction to determine the visibil- 
ity of the sources to be monitored in a microlensing 
survey. The flux density received from a source at £ 
is 

(A10) 



T = 



3 -0.921A /8 C 



The fraction of sources with L 1 > L I0 in Eqs. (A9) 
has the form C\L^ + C2. If T a is the minimum flux 
density that can be reliably monitored in a planet 
search from a source at D os , the cumulative luminos- 
ity function takes the form 



= a^,C,Ars) (All) 
= a (4nD 2 s F <; 2 e - 921A ^f + C 2 . 



Eq. (All) represents the fraction of sources that are 
visible at distance 

Although beyond current practice (P. Sackett, pri- 
vate communication, 2000), we shall assume a source 
with / = 21 can be monitored successfully from the 
ground in planet searches. This corresponds to the 
minimum usable flux density of J- — > J- 21 = 6.21 x 
10 _15 ergs cm _2 sec _1 . For the / band LF that we have 
adopted, about 1.16% of the stars have / 8 < 17 and 
about 4.18% have / 8 < 18.5. These are the points 
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in the LF where the slope changes, so the values of 
Ci, C2, f3 must change for F in Eq. (All) when 
F = 0.0116 and 0.0418. Of course F = 1 for suf- 
ficiently small values of £ where all of the stars are 
visible, and F = for sufficiently large values of £ 
where no stars are visible, where "visible" here means 
/ < 21. For example, substitution of the values of the 
constants from Eqs. ([Ai3) along with the value of T 21 



into Eq. (All) yields 



.000259, 



F{T>T 21 ) = t(F ,t,A Ia ) 

= 1, C < 0.0892, 

= 0.240 (C 2 e°- 921 ^^)" a326 -n,()7l!. 
0.0892 < C < 0.9075, 

= 0.300 ( C 2 e 0-921A„C) "0-827 

0.9075 < C < 1.2025, 
= 0..482 (^0.921^-1.01° _ .0000275, 

1.2025 < C < 2.7677, 
= 0, 2.7677 < C (A12) 

as the fraction of sources visible at distance ( for 
A Ig = 3. For A la = 1 the expressions are identi- 
cal to those in Eqs. (A12) except the breaks oc- 
cur at C = 0.0965, 1.5536, 2.2477, and 6.4602. In 
this latter case the most luminous star in the LF 
could be seen far beyond the opposite edge of the 
Galaxy even if the absorption were to persist all along 
the line of sight. For A I8 — 5 the breaks are at 
C = 0.0833, 0.6739, 0.8648 and 1.8384. Fig. [jj shows 
the fraction of visible sources (/ < 21) as a function 
of £ for three values of A IS . 



B. Amplification of the source 

We adopt the analysis of Witt (1990) to deter- 
mine the amplification of a point source star (source) 
by the star-planet binary lens as a function of the 
relative angular positions of the source, lensing star 
of mass M (hereinafter, the lens) and the planet 
to. However, since M 3> to, we choose the opti- 
cal axis to pass through the lens and locate source 
and planet relative to this axis. All displacements 
of the projected positions on the lens plane are nor- 
malized by the Einstein ring radius of M, R E = 
y/AGMD 0L (D 08 - D OL )/(c 2 D os ), where G is the 
gravitational constant, c is the velocity of light and 
D OL and D os are the observer-lens and observer- 
source distances respectively. (R E is the radius of 



the axially symmetric ring image of a source that is 
directly behind the lens.) The lens equation is then 
written 

X = w- — + — , (Bl) 



w* p* — w* ' 
where A = ^+irj and w = x+iy are the displacements 
of the source and image from the optical axis respec- 
tively, p — p x + i p y is the displacement of the planet 
and e = q. The (*) indicates complex conjugation. 
A, w, p may be interpreted as either the projected dis- 
placements in the lens plane normalized by R E or the 
angular displacements normalized by the angle sub- 
tended by R E . The axes (£,77) and (x,y) are aligned. 



If we multiply the complex conjugate of Eq. (Bl) 
by the product w(p—w) and substitute w as a function 
of w* from Eq. (Bl), the following polynomial results: 



P(w*) = a$w* +04UI* +a^w* +d2W* +aiw*+ao = 0, 

(B2) 

where 



a>5 

as 

a-2 

ai 
«o 



A(A - p) 

A(A(-A* - 2p*) + p\* + 2pp* + l + e)-ep 
A*(A(2p*A - 2(1 + e) - 2pp*) + p(l + e)) 
+X(Xp* 2 - pp* 2 - 2p*) + epp* 
X*(X(-p* 2 X + p*(pp* + 2e + 4)) - (2 + e)pp* 
-(l + e) 2 ) + Xp* 2 (l-e) + p*e(l + e) 
X*{-2p* 2 X + p*(pp* + 2(1 + e))) - ep* 2 



-AV 



(B3) 



The roots of Eq. (B2) are the complex conjugates 



of the positions of the images, where the number of 
valid images is either 3 or 5 with validity determined 
by substitution of the corresponding root into Eq. 



Eq. (Bl) represents a transformation of the coor- 
dinates of an image (x, y) to the coordinates of the 
source (£, rf). The amplification of the source repre- 
sented by a particular image is the ratio of an element 
of area at the image position in the (x, y) system, 
as transformed from the source position in the (£, r/) 
system, to the original element in the (£, 77) system. 
Hence, the total amplification of the source due to the 
star-planet binary lens is 



.4 



where (Witt, 1990) 



X]|det J, 



(B4) 



<9£ dr\ 
dx dy 



det(J) = ^ ^ 



dy dx 
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dX 
dw 



d\ ( dX V 



dw* \ 



= 1 



(p* — to*) 5 



(B5) 



is the determinant of the Jacobian of the transfor- 
mation represented by Eq. (Bl), and Ji = J(wi) is 
the Jacobian evaluated at the coordinates of the ith 
image. The sum is over either three or five images. 
For a single lens without the planet the amplification 
simplifies to 

u 2 + 2 

A = . Z (B6) 
u\/u z + 4 

where u is the angular separation of the lens and 
source normalized by the angle subtended by the Ein- 
stein ring radius in the lens plane. 

The curves defined in w space by det J — are 
called the critical curves where the amplification is 
formally infinite. This singularity is removed when 
the finite size of the source is considered. The cor- 
responding closed curves in A space found from Eq. 
(Bl) are the caustics, such that when the source is on 
a caustic curve an image is on a critical curve. For the 
single star lens, the caustic reduces to a single point at 
the lens position, and the corresponding critical curve 
is the Einstein ring. If the source never crosses a caus- 
tic during the event there are always just three im- 
ages and therefore three terms on the right-hand-side 
of Eq. (B4). At the first caustic crossing, an image 



is created on the critical curve and separates into two 
images for a total of 5 images that persist as long as 
the source is inside the closed caustic curve. On the 
second caustic crossing, with the source now outside 
the closed caustic, one of images created at the first 
crossing merges with one of the original three images 
and disappears — leaving again only three images. 

To evaluate contours of equal amplification for a 
given position of the source, the amplification is eval- 
uated for a 401 x 401 grid of planet positions centered 
on the lens and extending ±2.5i? E in both orthogonal 
directions. A 401 x 401 array of the fractional change 
in the amplification due to the planet, A/Aq — 1, is 
saved, where A is the amplification for a particular 
position of the planet and Ao is the unperturbed am- 
plification for the given source position. Contours of 
constant fractional change in the amplification (the 
perturbation due to the planet) are obtained with any 
contouring program for arbitrary fractional changes. 
In Fig. [2| the contours shown correspond to the as- 
sumed minimal perturbation that can be detected ac- 



cording to the S/N at that position on the light curve. 
However, these contours are not used to trace their ex- 
tremes along the line containing the source and lens 
during the course of the event. 

It is more efficient to determine these extremes di- 
rectly by moving the planet along the axis containing 
the source and lens to determine the two points on op- 
posite sides of each unperturbed image position where 
the fractional change in the amplification is 2 /(S/N) 
for the major image outside R E and —2/(S/N) for the 
minor image inside R E . The lens-source line is chosen 
to be the real axis so that A = A* and p = p* . The 



coefficients in Eq. (B3) are then real and analytic in 
p, which is now the linear position of the planet along 
the real axis. We form the function 



f(p) 



(B7) 



where A = Aq[1 ± 2 /(S/N)] defines the extremes of 
the planet positions along the real axis on cither side 
of the unperturbed image positions inside of which the 
fractional perturbation would exceed 2/ (S/N) and al- 
low detection of the planet. These extremes are de- 
fined by the zeros of f(p). A trial value of p somewhat 
removed from the unperturbed position of a real im- 
age is used to start a "safe" Newton's method solution 
for a zero of f(p) (Press, et al. 1986). All the nec- 
essary derivatives can be found from Eqs. (B2) and 



(B3), where p* = p in the latter, and the realness of 
the coefficients is exploited. A step in the procedure 
is thus for a trial position of p, all of the valid roots of 
the polynomial form of the lens equation (Eq. (B2)) 



are used to form f(p) and df(p)/dp for use in New- 
ton's method to determine the next approximation to 
the root, where the procedure is iterated until a tol- 
erable precision is obtained. The two pairs of roots 
are determined for each position of the source as it is 
stepped through the event, where S/N is a function 
of the source position. 



C. Probability distribution of projected star- 
planet separation 

If dn L (M, D OL )/dM is the number density of stars 
(lenses) per unit mass interval near mass M and near 
distance D OL , then (dn L (M, D OL )/dM)Dl L dD OL dMAfl 
is the number of stars in mass range dM about M and 
in volume D^ L dD OL Afi, where Af2 is a representative 
solid angle in the telescope field of view. The fraction 
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of the lenses in the volume of those along the LOS to 
the source is just this quantity divided by its integral 
along the LOS and over the mass function, where we 
assume the MF to be spatially invariant. Hence, 



o Jm 



F(x p ,a,() = 



F(x p ,a,(,z,M)Q 



a 2 C 



KMz(C-z) p 



z2 dn L (M,z) dMdz 



dM 



,dn L (M, z) 



dM 



dMdz 



(CI) 

lo Jm. 

is the average of F(x p , a, £, z, M) over the lens distri- 
bution from observer to source and over the MF. We 
have normalized the distances with D$ as above and 
introduced the step function ©(a;) to exclude those 
regions of z, M space for which the integrand in the 
numerator is undefined. If x 2 > a 2 (/(KMz(( — z)), 
r > a in Eq. (||). For fixed x p and a, Eq. (CI) 
represents the average of F(x p , a, £, z, M) over those 
lenses distributed along the LOS to distance £ for 
which R E < a/xp. 

For the density distribution in the Galaxy we adopt 
the double exponential model for the disk distribution 
of mass and a triaxial bar distribution for the bulge. 
The disk model is that of Bahcall and Soneira (1980), 



Pm 



M„ 



dp 
dM' 



dM 



Pmo exp 



( -\A r 
\300pc s a 



(C2) 

where z' is the coordinate perpendicular to the plane 
of the galaxy, r is the radial coordinate in the plane 
of the galaxy, and s d — 2.7 kpc is the scale length 
in the radial direction, which is chosen less than the 
3.5 kpc used by Bahcall and Soneira (Zhao, Spergel 
& Rich 1995; Kent, Dame & Fasio 1995), and where 
p M0 is chosen such that p M — O.O5Af /pc 3 at z' = 
and r = 8 kpc. Like Zhao, Rich and Spergel (1996), 
we choose the triaxial bulge model of Zhao (1996) 
but keep the nucleus in a truncated form and do not 
terminate the bulge at 3.3 kpc: 



Pm = Pa 



where 



exp 



-1.85 



exp(— s a ) 



(C3) 



fi2 = g a (s +y ) + z . f ^-j— j >q 

Of 

0.3136^ + z' 2 /— x^. neo , 

= ^ if \Jx 2 + y 2 < 0.56 kp<C4) 



with q a = 0.6, ct x = 1.49 kpc, a y = 0.58 kpc, and 
<r a = 0.40 kpc. The coefficient p Q determines the mass 
of the bulge. The nucleus is truncated to avoid the 
singularity at the origin in the Zhao bulge model. The 
observations are consistent with the long axis of the 
bulge inclined about 13° to 20° relative to the LOS to 
the Galactic center with the near side of the bar lying 
in the second quadrant. We shall choose 9 = —13°, 
but note that this angle is very uncertain. 

The origin of coordinates for both disk and bulge 
distributions is transferred from the Galactic center to 
the position of the observer with the following trans- 
formation: 

x — cos 9 — z cos b cos (£—6), 
y = — sin# — z cos 6 sin (I — 9), 
z' = zsinfr, (C5) 

for a lens in the bulge where z is the distance along 
the line of sight normalized by D s with xyz' being 
galactocentric coordinates along the bulge principal 
axes normalized by Ds and 9 the inclination of the 
bulge axis to the line between the observer and the 
Galactic center measured in the direction of increas- 
ing I. For the disk distribution, the transformation 



1 + z 2 cos 2 b ~ 2z cos b cos £, 



z' = z sin b, 



(C6) 



For a dynamically constrained bulge mass of 2.2 x 
10 10 M o (Zhao, 1996), p in Eq. (|c|) is 2.5M pc- 3 . 
The coefficient p M0 — 0.05e 8 / Sd M Q pc -3 (s d in kpc) 



in Eq. (|C2j) to satisfy the constraint on the local mass 
density. The conversion of these mass distributions to 
number density distributions for use in dn L /dM de- 
pends on the assumed mass function. We shall adopt 
the Holtzman et al. (1998) MF for both disk and 
bulge stars given by 



4>(M) = 0.7795Ar 
= 0.5081M" 



0.08 < M < 0.7 
0.7 < M < 2.0.(C7) 



The MF in Eq. (jCTj) is defined so that M<f>(M)dM 
represents the mass of stars/pc in mass range dM 
about M normalized such that the total mass density 
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is lM pc" 3 . A more recent MF (Zocalli et al. 2000) 
has an index of 1.3 instead of 1 in the M star region. 
However, we use the Holtzman et al. MF and point 
out the effects of increasing the index for small mass 
stars in Fig. [l7]. We truncate the MF at 2M be- 
cause of the small number of stars of higher mass. We 
could extend the MF into the brown dwarf region, but 



the microlensing event time scale frequency distribu- 



tion does not seem compatible with a large number of 
brown dwarfs — at least in the Galactic bulge (Peale, 
1999). 

The number of stars/M , or equivalently, stars/pc 3 
for this assumed MF is found from integration of Eq. 
(C7) to be 2.1562 stars/M , so the number density 



Albrow MD, An J, Beaulicu J-P, Caldwell JAR, Dc- 
poy DL, Dopminik M, Gaudi BS, Gould A, Grccnhill 
J, Hill K, Kane S, Martin R, Menzies J, Naber RM, 
Pel J-W, Pogge RW, Pollard KR, Sackett PD, Sahu 
KC, VErmaak P, Vreeswijk PM, Watson R, Williams 
A. (2000b) Limits on the abundance of galactic plan- 
ets from five years of PLANET observations 
ph/000807S| . Astrophys. J. In press. 



ast ro- 



of sta rs c orresp ond ing to the mass densities given by 
Eqs. (C2) and (C3) are found by multiplying each by 
this factor . The expression for dn . ldM is the sum of 



thcl disk and bulge contributions at a given point , so 



Alcock C, Allsman RA, Alves DR, Axelrof TS, Becker 
AC, Bennett DP, Cook KH, Drake AJ, Freeman KC, 
Geha AM, Griest K, Lehner MJ, Marshall, SL, Min- 
niti D, Nelson CA, Peterson BA, Popowski P, Pratt 
MR, Quinn PJ, Stubbs CW, Sutherland W, Tro- 
maney AB, Vandehei T, Welch DL. (2000) The Ma- 
cho project: Microlensing optical depth towards the 
galactic bulge from difference image analysis, astro-ph 
000251C . Astrophys. J. In press. 



dn L (M,z)/dM = 

0.7795[0.1078e A + 5.3903(e"^ + s^^e'^/M, 
0.08 < M < 0.7 

= 0.5081[0.1078e A + 5.3903( e -^ + s^^e-^^/M 2 - 2 , 

0.7 < M < 2.0 (C8) 



where A is the exponent given in Eq. (C2) +S/sd but 
with the change of coordinates given in Eqs. (C6) 
and (|C5|) in all of A, s a , andsfc. As we shall be deal- 



ing in fractional distributions, only the ratio of the 
coefficients of the disk and bulge distributions will 
have any effect on the probabilities. The integrals in 
Eq. ( pip are evaluated with a Monte Carlo technique 
(Press et al. 1986) as a function of x p for given a. 
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